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SUMMARY 


The static transmission error is widely recognized to be the principal source of 
vibratory excitation arising from meshing gear pairs. In this report, a generalized 
three-component transmission error of meshing spiral bevel gears is defined. Equa- 
tions and computational procedures are developed for predicting the generalized 
transmission error components from their underlying causes. 

First, it is shown that the traditional one-component transmission error is 
incapable of describing the vibratory excitation of spiral bevel gears arising from 
nonconjugate meshing action without making arbitrary, generally unjustifiable assump- 
tions. A three-component generalized transmission error definition is introduced to 
rectify this deficiency. Contact between meshing teeth is assumed to take place in a 
plane surface containing the pitch element. Gears with spherical involute teeth 
exactly satisfy this assumption, and gears manufactured on the crown rack basis 
closely satisfy this assumption. 

Equations are derived yielding the three components of the generalized trans- 
mission error in terms of deviations of tooth surfaces from equispaced perfect 
spherical involute surfaces and tooth/gearbody elastic deformations arising from the 
three components of the generalized force transmitted by the gears. It is explained 
how to incorporate these equations into the equations of motion of a gear system. 

I 

Equations are derived for the three components of the generalized force trans- 
mitted by the gears which are valid whenever inertial effects of the meshing gears 
and their supports are negligible. Bearing offsets from the positions occupied by 
| the shaft centerlines of perfect spherical involute bevel gears and bearing/bearing 

i support flexibilities enter into the computation of these forces. When these forces 

are used to predict the transmission error, one obtains the generalized static 
transmission error. 

For given forces transmitted by the gears, the three components of the general- 
ized transmission error are shown to be stationary with respect to small independent 
variations in the positions of the endpoints of the lines of tooth contact about 

| their true values. One of these transmission error components is shown to take on a 

minimum value with respect to these variations. A computational procedure based on 
this extremum principle is described which is suitable for computing the zone of 
tooth contact in the plane of contact, the three components of the generalized static 
transmission error, and their spectra. 



INTRODUCTION 


For some time, lack of conjugate action (refs. 1,2) has been recognized to be 
the principal source of vibratory excitation arising from the meshing action of 
parallel-axis gear pairs (refs. 3-11). This lack of conjugate action is approxi- 
mately characterized by the static transmission error (refs. 3,4), which is loosely 
defined (refs. 7,9) as the deviation from linearity of the angular position of a gear 
measured as a function of the angular position of the gear it meshes with when the 
gear pair is transmitting torque at low enough speed so that inertial effects are 
negligible. A more precise definition of the static transmission error is given by 
equation (1) of reference 9 and equation (2.1) of this report. 

The static transmission error can be combined with the equations of motion of a 
gear system to predict responses of system elements of interest using the approximate 
method described on p. 1410 of reference 7 and pp. 13-18 of reference 12. An exact 
method for combining the transmission error with the quations of motion is described 
by equations (3.44) and (3.45) and the accompanying text of this report. 

When characterizing the vibratory excitation of parallel-axis gear systems by 
the transmission error, there is an implied assumption that the shafts in a meshing 
gear pair remain parallel with one another even though they may be moving in relation 
to each other (refs. 7,9,12). Such motion is necessary for vibratory energy to be 
transmitted from the shafts into the bearings and gear casing - which, in turn, 
vibrates and radiates sound. However, in the case of nonparallel-axis gears, and 
spiral bevel gears in particular, it is shown in the present report in the text 
accompanying equations (2.3) and (2.4) that there would appear to be no physically 
justifiable assumption comparable to the assumption of the shafts of parallel-axis 
gears remaining parallel during their motions. Thus, in the present report, the 
traditional concept of the transmission error is generalized to a three-component 
transmission error. The three components completely characterize the vibratory 
excitation arising from arbitrary relative motions of the two gears of a meshing 
pair - provided that the surface of tooth contact is reasonably well approximated by 
a plane surface passing through the point occupied by the pitch cone apexes of the 
two meshing gears. 


2 



The generalization of the traditional one-component transmission error to the 
three-component generalized transmission error very substantially complicates trans- 
mission error analysis and its usage in the equations of motion of gear systems. 
Nevertheless, this generalization is necessary to avoid completely arbitrary assump- 
tions characterizing the angular motions between nonparallel shafts. 

However, the three-component generalized transmission error described herein 
retains another assumption employed in the traditional transmission error char- 
acterization; namely, that all tooth contact takes place in a plane surface. In the 
case of parallel-axis gears, this plane is the base plane (refs. 7-9), sometimes 
referred to as the plane of action. In the case of straight and spiral bevel gears, 
this surface is the base plane associated with spherical involute tooth generation 
(p. 199 of ref. 13) which is discussed in detail in the present report. The base 
plane also may be regarded as the plane containing the pitch element that is tangent 
to the true surface of tooth contact (pp. 1-4 and 1-16 of ref. 14). The octoid path 
of contact of gears manufactured on the crown rack basis (pp. 494-496 of ref. 15) is 
very nearly linear where actual tooth contact takes place (p. 322 of ref. 1; p. 1-15 
of ref. 14). Thus, it is felt that the very substantial simplification achieved by 
assuming all tooth contact to take place in a plane surface containing the pitch 
element easily justifies this basic assumption. Even utilizing this assumption, the 
analysis required to understand and predict the generalized transmission error is 
very much more complicated than that required to understand and predict the tradi- 
tional transmission error of parallel-axis gears. 

The goals of the work described in this report have been to carry out an 
analysis of the generalized transmission error arising from spiral bevel gear pairs 
in sufficient detail (1) to gain a physical understanding of the factors that govern 
its magnitude and characteristics, and (2) to predict its magnitude and character- 
istics from these factors it depends on. It is shown that these physical factors 
include: 

(1) deviations of tooth running surfaces from perfect conjugate surfaces 

(2) tooth/gearbody stiffnesses 

(3) deviations of gearshaft bearing centerline positions from the shaft centerline 
positions of perfect conjugate gears 
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(4) bearing/bearing support flexibilities, and 

(5) the shaft input torque. 

Appropriate descriptions of all of the above quantities are required to predict 
the generalized static transmission error of spiral bevel gears. Thus, all of the 
above quantities play important roles in controlling the magnitude and characteris- 
tics of the vibratory excitation arising from spiral bevel gear pairs. 

In addition to the above-listed factors, the contact ratios play an important 
role in smoothing the lack of conjugate action contribution of individual teeth 
(refs. 7-9). Large contact ratios are helpful in reducing the transmission error. 
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THE GENERALIZED TRANSMISSION ERROR OF SPIRAL BEVEL GEARS 


It is widely recognized (refs. 3-11) that the static transmission error (refs. 
3,4) describes the principal source of vibratory excitation arising from meshing 
gear pairs. The static transmission error characterizes the deviation from 
| conjugate action (refs. 13,16) of one gear with respect to its mating gear. The 
i angular static transmission error can be defined loosely (refs. 7,9) as "the 
deviation 60 from linearity of the angular position 0 of a gear measured as a 
function of the angular position of the gear it meshes with when the gear pair is 
| transmitting a constant torque at low enough speed so that inertial effects are 
| negligible." In gearing vibration analyses, it generally is more convenient to use 
a lineal equivalent 5 = R^60 to the angular static transmission error, where is 
the base circle radius of the gear whose angular transmission error is 60. The 
lineal static transmission error describes the deviation from conjugate action as a 
linear error in the plane of contact. 

A slightly different definition of the lineal static transmission error useful 
for characterizing the vibratory excitation arising from a meshing gear pair of 
I nominal involute design can be given (ref. 9, p. 57) as 


4 .( 1 ) ..( 1 ) .( 2 ) „„( 2 ) 


5 =R b 


60 


- R. 


60 


( 2 . 1 ) 


where R^^ and R^^ are the base circle radii of the two gears (1) and (2) of the 
pair and 60^^ and 60^^ are the angular deviations of gears (1) and (2) from the 
positions of their rigid perfect involute counterparts determined under the condi- 
tion that lateral displacements v^^ and v^) of the two gear shafts in the 
direction parallel to the plane of contact are zero (Fig. 1). The sign convention 
implied by equation (1)* and Fig. 1 is such that the lineal static transmission 
error 5 is positive when the deviations from perfect involute teeth are equivalent 

j to removal of material from the running surfaces of the teeth. 

I 

I 

j In a gearing system in operation, the transmission error from each gear pair in 

j the system arising from static and dynamic effects must be compensated for by rota- 
j tional increments and torsional deformations of the shafting and/or lateral 


1 *Reference to equations in the same section as the equation citation will be made 
without the section number in the citation. 

! 
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translational displacements v^^ and v^) of the gear shafts in the direction 
parallel to the plane of contact (Fig. 1). This observation can be expressed (ref. 
7, p. 1410) by the relation 


g(l) 5 » 0 ( 1 ) + V (D 



6'e 


(2) _ v (2) 


= C 


( 2 . 2 ) 


where the angular deviations 6 * 0 ^^ and 6 ' 0 ^) encountered under running conditions 
generally are different from the angular deviations 60^^ and 50 respectively, 
used in the definition ( 1 ) of the static transmission error. 

In utilizing the static transmission error 5 as the vibratory excitation in the 
equations of motion of a gearing system, an equation of the form of equation ( 2 ) is 
written for each pair of parallel-axis gears in the system. Since the lateral shaft 
displacements v^) and v^) are required to predict the rate of vibratory energy 
flow into the gear supporting structure, accurate prediction of these displacements 
is very important in noise analyses. However, in equation (2) there is an implicit 
assumption that the shafts of gears ( 1 ) and ( 2 ) remain parallel when laterally 
displaced - i.e., that each displacement v^^ and v^) is the same at the bearing 
locations at the two ends of the shaft. The simplification of the analyses that 
results from this assumption undoubtedly justifies its use in most cases involving 
parallel-axis gears. 

In the case of bevel gears, however, there does not seem to be a simplifying 
assumption of comparable validity concerning lateral shaft displacements that would 
lead to a single relation analogous to equation (2) for each gear pair. To see 
this, we note first that, for either gear ( 1 ) or gear ( 2 ), the transmission error ? 
in equation ( 2 ) can be compensated for by the term r£* ^ 6 ' 0 ^*^ arising from the 
angular deviation 6 ' 0 ^*^ or by the term v^* ^ arising from the lateral shaft 
displacement of gear (•). For parallel-axis gears, the base circle radius is a 

constant that is independent of axial location on the gear. In the case of bevel 
gears of nominal (spherical) involute design (ref. 13, p. 199; ref. 16, p. 128), 
however, the base cylinder of a parallel-axis gear is replaced by a base cone whose 
base circle radius varies in direct proportion to axial distance from the base cone 
apex. Thus, a single equation of the form of equation (2) could be made to apply to 
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bevel gears only if one were to assume that the quantity £-v^^ + varies 

axially in direct proportion to distance x' from the base cone apex - i.e., that 

? - v (1) + v (2) = ax', (2.3) 

where a is a constant. If the transmission error t, is assumed to be independent of 
axial location x', this assumption would, in turn, require that shaft lateral dis- 
placements v^D and v^) vary with axial location x' as 

v (1) - v^ = ? - ax', (2.4) 

where 5 and a are independent of x'. It would be difficult to justify such an 
assumption on physical grounds. 

In the present section, the need for such an assumption is removed by defining 
the generalized transmission error for spiral bevel gears of nominal spherical 
involute design. The generalization of the concept of the transmission error 
contained herein removes the need for all assumptions in gearing vibration analyses 
that would constrain the displacements of gear shafts from the positions occupied by 
the shafts of their rigid perfect involute counterparts. 

Kinematics of Involute Spiral Bevel Gear Action 

Geometric Description of Tooth Contact. Straight-tooth and spiral bevel gears 
transmit motion between nonparallel intersecting shafts. The motion transmitted by 
perfect bevel gear pairs is equivalent to that transmitted by perfect right circular 
cones with a common apex in pure rolling contact with no slipping (ref. 13, pp. 125- 
128). The ratio of the instantaneous angular velocities of the two cones about 
their axes is exactly constant. Such cones in rolling contact are referred to as 
pitch cones. 

The teeth of a pair of involute bevel gears may be imagined to have been 
generated utilizing the base cones of the gears. The base cones are smaller than 
and concentric with the pitch cones. The base and pitch cones of both meshing gears 
have the same common apex A (Fig. 2). Let us imagine that the base cone of the 
upper gear in Fig. 2(a) is covered by a thin inextensible membrane which is un- 
wrapped without slipping from the upper gear and rewrapped without slack onto the 
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CONE 

(a) (b) 


FIG. 2. SPHERICAL BASIS OF BEVEL GEAR ACTION (ADAPTED FROM FIG. 5.18 OF REF. 17 
WITH PERMISSION). 
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base cone of the lower gear as the gears rotate. A point on this membrane traces 
out a spherical involute curve (ref. 13, pp. 197-199) as it is unwrapped from the 
base cone of the upper gear. As the membrane is unwrapped, such a point may be 
imagined to trace out the running suface of the tooth t as is indicated in Fig. 2(a). 
This process of generating a spherical involute is entirely analogous to the process 
of generating the involute of a circle (ref. 13, p. 163-167), except that the 
involute of a circle is generated by unwrapping a membrane from a right circular 
cylinder, whereas a spherical involute is generated by unwrapping a membrane from a 
right circular cone. 

As the membrane is unwrapped from the base cone of the upper gear shown in Fig. 
2(a) and rewrapped onto the base cone of the lower gear, a point on the membrane 
traces out a trajectory that lies on a spherical surface - as is indicated in the 
figure. Conversely, when the membrane is unwrapped from the cylinder in the genera- 
tion of an involute of a circle, a point on the membrane travels in a plane that is 
perpendicular to the axis of the cylinder. This plane is called the transverse 
plane. Thus, the spherical surface with center at the base and pitch cone apex 
illustrated in Fig. 2 is the involute bevel gear counterpart to the transverse plane 
of involute parallel-axis gears (ref. 17, pp. 56-58). We shall call this spherical 
surface the transverse sphere . 

Let us now imagine several concentric spherical sections of the type shown in 
Fig. 2, each located at a different radius from the common apex of the base and pitch 
cones. On each concentric sphere, the teeth of each of the two meshing gears are 
assumed to have been constructed as spherical involutes as described above. The loci 
of points connecting the spherical involute surface of a generic tooth t from one 
sphere to the next determine the shape of the running surface of that tooth. In the 
case of straight-tooth bevel gears, these loci of points are straight lines normal to 
the spheres which originate at the apex of the base and pitch cones. Straight-tooth 
bevel gears are the bevel gear counterparts to spur gears. In the case of spiral 
bevel gears, the loci of points connecting the surface of tooth t from one sphere to 
the next are curved or spiral lines which also originate at the apex of the cones, 
but which are not generally normal to the spheres. 

Any spiral bevel gear may be imagined to have been generated from a straight- 
tooth bevel gear by cutting the straight-tooth gear into an (infinite) number of 
transverse spherical sections and rotating the sections relative to one another about 
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the axis of rotation of the gear. This operation produces the stepped bevel gear 
counterpart (ref. 13, p. 202) to the conceptual generation of a parallel-axis helical 
gear from a spur gear, which is carried out by cutting the spur gear into an 
(infinite) number of sections normal to the gear axis and rotating the sections 
reltive to one another (ref. 13, pp. 190,191). For any shape of tooth spiral, the 
kinematic action of an involute spiral bevel gear is identical in any transverse 
spherical section to that of an involute straight-tooth bevel gear. 

The taut inextensible membrane connecting the base cones of the two meshing 
bevel gears (Fig. 2) is tangent to the conical surfaces of both base cones. The 
locus of points of tangency between each base cone and the membrane is a straight 
line. From the fact that these two straight lines of tangency intersect at the 
common apex of the two base cones it follows that the portion of the taut membrane 
lying between the two base cones is a plane surface tangent to both base cones. 
Furthermore, from the spherical involute construction of the tooth surfaces it 
follows that all tooth contact takes place in this plane of tangency for straight- 
tooth and spiral bevel involute gears. This plane of tangency is called the base 
plane (ref. 14, pp. 1-16). 

At this juncture, it is convenient to replace the taut membrane by a rigid 
circular disk of infinitesimal thickness positioned in the plane of contact with 
center at the common apex of the base cones . The disk lies in the base plane and is 
tangent to both base cones (Fig. 3). As the two gears rotate, the base disk rotates 
about its center without slipping on either of the two base cones. The motion of the 
portion of the disk that momentarily lies between the lines of tangency to the base 
cones is identical to the motion of the taut inextensible membrane that it replaces. 

Since all tooth contact takes place in the base plane, all instantaneous lines 
of tooth contact occur in the plane of the base disk. It follows from the spherical 
involute construction of the tooth surfaces that the circumferential distances 
separating the lines of tooth contact in the base plane do not change as the gears 
rotate. We may imagine that these lines of contact are drawn on the base disk. As 
the gears and base disk rotate about their respective axes, the true lines of tooth 
contact and the lines of contact drawn on the base disk exactly coincide in the 
region of the base plane where tooth contact occurs. 
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BASE DISK 
PLANE OF CONTACT 

BASE 

CONE 

PITCH 

CONE 


ADDENDUM 

CONE 


FIG* 3. LOWER FIGURE SHOWS BASE DISK TANGENT TO BASE CONES OF TWO MESHING PERFECT 
SPHERICAL INVOLUTE SPIRAL BEVEL GEARS. UPPER FIGURE SHOWS ZONE OF CONTACT 
IN BASE PLANE AND LINES OF CONTACT ON BASE DISK. (LOWER FIGURE ADAPTED FROM 
FIG. 5.18 OF REF. 17 WITH PERMISSION.) 
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Figure 3 also shows the addendum cones of the two gears which share a common 
apex with the base and pitch cones. Actual tooth contact takes place in the sector 
of the base plane located between the intersections of the two addendum cones with 
the base plane as is shown in Fig. 3. The two intersections delineating the zone of 
contact are straight radial lines that intersect at the common apex of the addendum, 
base, and pitch cones (Fig. 3). 

The bevel gears are located between outer and inner spherical surfaces of the 
type illustrated in Figs. 2 and 3. If the spherical surfaces shown in Figs. 2 and 3 
are taken to be the outer surface enclosing the two gears, then the zone of tooth 
contact in the base plane is the portion of the sector between the outer and inner 
radii drawn on the base disk in Fig. 3. 

The lines of tooth contact in Fig. 3 are drawn solid within the zone of contact 
and dashed outside the zone of contact. As the gears and base disk rotate about 
their respective axes, the lines of tooth contact affixed to the base disk pass 
through the fixed zone of tooth contact as can readily be imagined with the aid of 
Fig. 3. This geometric model of spiral bevel gear action is an exact representation 
of the kinematic behavior of the lines of contact of perfect involute spiral bevel 
i gears. 


Tooth Contact Coordinates in the Base Plane. Concentric circles located at 
' radii r 0 + (F/2) and r 0 “ (F/2) on the base disk in Fig. 3 define the outer and inner 
I radii of the zone of tooth contact in the base plane. Thus, the radial span of the 
zone of tooth contact in the base plane is F and the radial midpoint of this span is 
! r Q . Let y denote a radial coordinate on the base disk in Fig. 3 defined by 



|l where r is the usual radial coordinate with origin at the center of the disk. From 
equation (5), it follows that the origin of y is located at r ■ r Q and that y is 
positive for values of r less than r Q . 

I 

i Lines of tooth contact on the base disk in Fig. 3 are designated by index values 

; j ■ 0,1,2,»»». As the gears and base disk rotate about their respective axes, the 
angular position of tooth pair j = 0 relative to the center of the zone of contact is 
designated by 0, which is measured in radians. When the angular position of the line 
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of contact of tooth pair j = 0 at radial location y = 0 exactly coincides with the 
centerline of the zone of contact on the base disk, we have 0 = 0 as is readily 
apparent from Fig. 3. 

Instead of dealing with the angular variable 0, we shall deal primarily with the 
lineal variable 

x = r 6, (2.6) 

o 

where r Q and 0 are defined above. Thus, x also determines the angular position of 
tooth pair j = 0 relative to the zone of contact centerline. Notice from Fig. 3 that 
x is the arc length in the base plane of the circular segment of radius r Q between 
the line of contact of tooth pair j = 0 and the centerline of the zone of contact. 

The tooth spacing A on the base disk in Fig. 3 also is measured on the same circular 
arc of radius r Q » 

Consider the transverse spherical section that passes through the point P shown 
in the zone of contact on the base disk in Fig. 3. The radius of this spherical 
section is uniquely determined by the coordinate y that describes the position of P 
on the line of tooth contact in the base plane. For the particular value of x = 
r Q 0 illustrated in Fig. 3, the point P determines a unique point located on the 
running surfaces of the two mating teeth. It follows that specification of a pair of 
coordinate values x,y determines a unique point of contact on the running surfaces of 
two mating teeth designated by index j. 

Let 0^^ and 0^^ denote, respectively, the rotational positions of gears (1) 

and (2) about their axes of rotation. Consider infinitesimal increments d0 , d0^\ 

( 2 ) 

and d0 v ' in the rotational angles of the base disk and the two meshing gears. From 
the fact that the base disk may be imagined to be rotating without slipping on the 
two base cones illustrated in Fig. 2, it follows that the infinitesimal lineal incre- 
ment dx is related to increments d0 , d0^\ and d0^) by 

dx = r^ d0 = RO) d9 (0 , *C2> d9 <2> , (2.7) 

where and are the radii of the base cones of gears (1) and (2), 

respectively, that join the base cone generatrices at distances r Q from the common 
apex as illustrated in Fig. 2(b). The integral of equation (7) yields 
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( 2 . 8 ) 


R< 2 V 2 > . 

Dividing equation (8) by r Q and using elementary trigonometric relations arising from 
Fig. 2(b) yields, further, 

(x/r Q ) * 9 * (siny^ 1 ^) 0^0 = (siny£ 2 ^) 8^ , (2.9) 

which relates the angular positions of the base disk and the two meshing gears using 
only the base cone angles yj^ and y£ 2 ^ shown in Fig. 2(b). 

Since yf^ or y^ 2 ^ can be chosen more or less arbitrarily, it follows with 

D D 

reference to the base disk in Fig. 3 that 2irr o /A need not necessarily be an integer. 
Thus, once 9 exceeds 2ir , the positions of the lines of tooth contact on the base disk 
do not, in general, repeat. 

The Generalized Transmission Error 

We now consider deviation of the kinematic behavior of real bevel gear pairs of 
nominal spherical involute design from the above-described ideal behavior of rigid 
perfect involute gear pairs. The effects considered are those arising from geometric 

I 

! deviations of the unloaded tooth surfaces from perfect spherical involute surfaces 
and from elastic deformations of the teeth and gear bodies caused by tooth loading. 

In the case of parallel axis spur and helical gears, the above effects have been 
characterized (refs. 3-11) by the transmission error, which is a one-component scalar 
representation. 

A more complete characterization of the above effects is provided by a three- 
component version of the transmission error, which we shall call the generalized 
transmission error. Let us consider the real gear counterparts to the perfect 
involute gears illustrated in Fig. 3. For any given value of x or 0 , there exists a 
lineal force distribution along each of the lines of contact within the zone of 
contact shown in Fig. 3 which represents the continuous distribution of forces 
transmitted by the pair of teeth in contact. If we neglect the effects of friction 
between meshing teeth, then all such distributed force vectors are normal to the 
tooth surfaces and from Involute construction therefore must lie in the plane of 
contact. The resultant total force transmitted by the mesh is the vector sum of all 


r 0 « r(1)q(1) 
O DO 
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forces transmitted by meshing tooth pairs; hence, the total transmitted resultant 
force vector also lies in the plane of contact. 

This resultant force vector can be described by three scalar components - 
namely, two orthogonal force components that lie in the plane of contact, say W t and 
W and the moment, say t, of the resultant total transmitted force taken about a 
specified point in the plane of contact. The directions of W fc and W r will be chosen 
perpendicular and parallel to the zone of contact centerline, respectively. Moment 
x will be taken about the point located at the mid-radial position y = 0 on the zone 
of contact centerline illustrated in Fig. 3. Force components W fc and W r will be 
taken as positive when they result in tooth pairs being brought together (rather than 
separated). If the total transmitted resultant force vector is taken to be positive 
for this same condition, then t will be taken as positive when the total resultant 
force vector passes through the zone of contact centerline at a negative value of y = 
r Q ~ r. This sign convention is illustrated in Fig. 4. 

The above-described distributed tooth forces give rise to elastic deformations 
of the teeth and bodies of the two meshing gears. These elastic deformations result 
in displacements of the gear bodies ''measured” at the gear shaft centerlines. 
Geometric deviations of the tooth running surfaces of the unloaded real gears from 
the spherical involute surfaces of their perfect involute counterparts give rise to 
additional displacements of the gear bodies at their shaft centerlines from the 
positions of their perfect involute counterparts. For a given value of x or 0 and 
given loading components W t , W r , and T , these displacements are characterized by the 
generalized transmission error which describes the composite change in position of 
the bodies of two meshing real gears at their shaft centerline locations from the 
positions of their rigid perfect involute counterparts. 

In the present report , the total deviation of the loaded tooth surfaces of the 
real gears from the surfaces of their rigid perfect involute counterparts is assumed 
to be sufficiently small so that the surface of actual tooth contact can be assumed 
to be the plane of contact of their rigid perfect involute counterparts illustrated 
in Fig. 3. Thus, in comparing the positions of the loaded real gears with the 
positions of their perfect involute counterparts, we are concerned only with those 
differential position components that lie ^n_ the plane of contact. It follows that 
the generalized transmission error can be characterized by three displacement 
components - namely, two translation components that lie in the plane of contact and 
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FIG. 4. SCALAR CHARACTERIZATION OF TOTAL FORCE ft TRANSMITTED BY THE GEAR MESH. THE 
FIGURE ILLUSTRATES INSTANTANEOUS LINES OF TOOTH CONTACT AND ZONE OF CONTACT 
IN PLANE OF CONTACT WHICH IS PLANE OF PAPER. FORCE COMPONENTS SHOWN ARE 
POSITIVE AND ARE ACTING ON THE GEAR TO THE RIGHT. MOMENT t IS POSITIVE WHEN 
ft PASSES THROUGH ZONE OF CONTACT CENTERLINE AT A NEGATIVE VALUE OF y AS 
SHOWN IN THE FIGURE. 
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one rotation component about a specified axis in the plane of contact* Since these 
components represent small displacements, the rotation component of the generalized 
transmission error can be represented by a vector (ref* 18, pp. 107,108) which is 
normal to the plane of contact. 

Each of these three displacement components of the generalized transmission 
error characterizes the composite displacement of both gears contributing to that 
component. The axis of rotation of the rotation component p of the generalized 
transmission error will be taken to be the point 0_ located on the zone of contact 
centerline at the value of y * 0 shown in Fig. 4. The rotation p represents the sum 
of the rotations about the axis 0_ of the two meshing gears from the positions of 
their perfect involute counterparts; p is taken to be positive when the portions of 
the teeth in the negative y region of the zone of contact come together more, rela- 
tive to their perfect involute counterparts, than the portions of the teeth in 
positive y region of the zone of contact. The directions of the two translation 
components g and g^ of the generalized transmission error are chosen to coincide, 
respectively, with the directions of the corresponding force components W t and W r 
illustrated in Fig. 4. Each component g and g^ represents the sum of the dis- 
placements of the two meshing gears from the positions of their perfect involute 
counterparts; g and g^ each is taken to be positive when it represents a displace- 
ment of the real gears "equivalent" to that caused by removal of material from the 
tooth surfaces of their perfect involute counterparts. Thus, positive values of 
W and x are associated with positive contributions of the elastic deformation 
components of g , g , and p, respectively. The sign convention of the generalized 
transmission error components is illustrated in Fig. 5. 

Expression for Generalized Transmission Error in Terms of 
Gear Displacements at Shaft Centerlines 

The base cone and plane of contact of a generic bevel gear are shown in Fig. 6. 
Orthogonal components xj # \ x^*\ and designate the three translation components 

of the shaft centerline position of gear (•) from the position of its rigid perfect 
involute counterpart. The origin of the x^ # \ x 2*^> x 3 # ^ coordinate system is 
located at the common apex of the base cones of the involute counterparts of the two 
meshing gears. The Xg*^ axis coincides with the axis of rotation of the involute 
counterpart of gear (•). The x^* ^ axis lies in the plane of contact. Positive 
directions of the Xj*\ X 2 *\ and x| # ^ axes are indicated by the arrows on the ends 
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FIG. 5. COMPONENTS OF THE GENERALIZED TRANSMISSION ERROR. MUTUALLY APPROACHING 
ARROWS REPRESENT GEAR BODIES COMING TOGETHER RELATIVE TO THEIR PERFECT 
INVOLUTE COUNTERPARTS. DIRECTIONS SHOWN REPRESENT POSITIVE VALUES OF 
Z t > C r , AND p. EACH COMPONENT REPRESENTS THE SUM OF THE DISPLACEMENTS OF 
THE TWO MESHING GEARS RELATIVE TO THEIR PERFECT INVOLUTE COUNTERPARTS; THESE 
DISPLACEMENTS ARE "MEASURED" AT THE GEAR SHAFT CENTERLINES. 
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FIG. 6. BASE CONE AND PLANE OF CONTACT FOR A GENERIC BEVEL GEAR PAIR. COORDINATES 
AND i = 1,2,3 MEASURE THE DIFFERENCE IN POSITION OF THE SHAFT 

CENTERLINE OF GEAR (• ) FROM THE POSITION OF ITS RIGID PERFECT INVOLUTE 
COUNTERPART. SIGN CONVENTIONS FOR THE x^ AND i " 1,2,3 ARE 

ASSOCIATED WITH THE GEARS SHOWN IN FIG. 3 AND DESCRIBED IN THE FIRST TWO 
ROWS OF TABLE 1. 


20 


of the coordinate axes. These designations completely define the x[*\ X2 *\ x 3 * ^ 
coordinate system with respect to the location of the base cone of the rigid perfect 
involute counterpart of gear (•). The base cone of the mating gear lies to the left 
of the base cone shown in Fig. 6. 

In addition to the three translation components x£*\ Xj*^, and x 3 ’\ three 
rotation components are required to describe arbitrary displacements of a real gear 
from the position of its perfect involute counterpart. Since such changes in 
position always are very small, the rotation components can be treated as vectors 
(ref. 18, pp. 107,108). We shall take the three rotation components to be the 
(small) rotations Pi*\ P2 *^> an< ^ P3*^ about the base cone axes Xj*^, X 2 *\ and x 3 *^, 
respectively. Since the x[*^ axis lies in the plane of contact, small rotations 
Pj*^ about this axis produce displacements at the plane of contact that are normal to 
that plane. Such displacements contribute nothing to the generalized transmission 
error, as indicated earlier. Thus, we shall be concerned only with the small rota- 
tions pg* ^ and p 3 *^ about the axes x|* ^ and x 3 *\ 

To obtain the contributions of Xj*\ x^* \ and Xg* ^ to the generalized trans- 
mission error, we require the contributions of these translations to the components 
and 5^ in the plane of contact shown in Fig. 5. These contributions can be 
obtained by transforming x[*\ and x 3 *^ into the rectangular coordinate system 

l z 1 , z 2 , Zg illustrated in Fig. 6. The plane containing the zj and Z3 axes is the 

plane of contact; the z, axis coincides with the zone of contact centerline shown in 

1 

Fig. 5. This transformation is carried out by two rotations - first, a rotation 
about the Xj* ^ axis through the base cone angle ^ yielding displacements in 
the y^ , y 2 , y 3 coordinate system followed by a rotation about the y 2 axis through the 
angle v^* ^ yielding the desired displacements in the Zj , z 2 , z 3 coordinate system. 

The resulting displacements along the z^ axis contribute to 5 and the resulting 
1 displacements along the z 3 axis provide the contribution of gear (. ) to £ r . 

Displacements along the z 2 axis are normal to the plane of contact and therefore do 
not contribute to the generalized transmission error. 

The only rotation components contributing to the generalized transmission error 
are those components about the z 2 axis, since (small) rotations about the z\ or the Z3 
axis provide displacement components at the plane of contact that are normal to that 
plane. Small rotations and p 3 *^ both contribute to the contribution of gear (. ) 

to the generalized transmission error rotation component y; in addition, P2* ^ and 
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Pg * ^ both contribute to 5 , since is the generalized transmission error transla- 
tion component in the direction that passes through the z 3 axis a distance r Q from 
the base cone apex (Figs. 5 and 6). 

The generalized transmission error components represent the sums of the 
contributions, as described above, from each of the two meshing gears. Let us 
designate the generalized transmission error by a column matrix or vector 

y(x) 

C (x) 
t 

C (x) 
r 

and the rotational and translational components of the shaft centerline of gear (•) 
measured from the position of its rigid perfect involute counterpart by the 
generalized displacement column matrix or vector 




q (,) (x) ^ 


( Pi°(x) 
P 2 * ^(x) 
P 3 * ^ (x) 


Xj* \x) 


.(•) 


(x) 


lx3*\x)J , 


( 2 . 11 ) 


(•) * (1) or (2), as appropriate. It is shown in the appendix that ?(x) can be 
obtained from the generalized displacements q^\x) and q^\x) of the two meshing 
gears by the matrix multiplications and summation 


*(.) - T (1, q (1) (,) ♦ T (2 Y 2 >00 , 


where \ (•) = (1) or (2), is defined by 


( 2 . 12 ) 
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> 


(2.13) 

where c and s are abbreviations (ref. 18, p. 106) for cosine and sine operations, 
respectively, on the angles ^ and v^*^. These angles and r Q are illustrated in 
Fig. 6. 


The terms in the bottom row in the right-hand side of equation (13) have two 
signs. Determination of the appropriate sign is given in Table 1. The gear pair 
shown in Fig. 3 is represented by the first line in Table 1. If gear (2) in Fig. 3 
becomes the driver rotating clockwise, the resulting configuration is represented by 
the second line in Table 1. The gear pair shown in Fig. 7 illustrates the other 
possible gearing configurations. Figure 8 displays the coordinate system used for 
each of the two gears shown in Fig. 7. The gear pair shown in Fig. 7 is represented 
by the third line in Table 1. If gear (2) in Fig. 7 becomes the driver rotating 
counterclockwise, the resulting configuration is represented by the fourth line in 
Table 1. Definitions of spiral directions are the same as those of screw threads 
(ref. 19, p. 241). Equations (10) through (13) and Table 1 can be regarded as the 
mathematical definition of the generalized transmission error for spiral bevel gear 
pairs . 
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Table 1 The upper or lower signs In the third row in equation (13) are chosen by the 
rules listed below. Each full line in the table represents one of the four possible 
configurations of spiral bevel gear pairs. Gear notions are observed fron large end 
of base cone. The first two lines in the table are associated with Figs. 3 and 6 and 
the 3rd and 4th lines are associated with Figs. 7 and 8. 


Upper Sign 

Driver clockwise, 
left-hand spiral 

Follower counterclockwise, 
left-hand spiral 

Driver counterclockwise, 
right-hand spiral 

Follower clockwise, 
right-hand spiral 


Lower Sign 

Follower counterclockwise, 
right-hand spiral 

Driver clockwise, 
right-hand spiral 

Follower clockwise, 
left-hand spiral 

Driver counterclockwise, 
left-hand spiral 
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ADDENDUM 

CONE 


ADDENDUM 

CONE 

PITCH 

CONE 


BASE 

CONE 


FIG. 7. ALTERNATIVE GEARING CONFIGURATION TO THAT SHOWN IN FIG. 3. CONFIGURATION 

SHOWN ABOVE IS ASSOCIATED WITH BASE CONE AND PLANE OF CONTACT CONFIGURATION 
SHOWN IN FIG. 8 AND DESCRIBED IN THE LAST TWO ROWS OF TABLE 1. (LOWER 
FIGURE ADAPTED FROM FIG. 5.18 OF REF. 17 WITH PERMISSION.) 


25 





FIG. 8. BASE CONE AND PLANE OF CONTACT FOR BEVEL GEARS SHOWN IN FIG. 7 AND DESCRIBED 
IN THE LAST TWO ROWS OF TABLE 1. COORDINATES x^' ) AND p^ #) , i = 1,2,3 
MEASURE THE DIFFERENCE IN POSITION OF THE SHAFT CENTERLINE OF GEAR (•) FROM 
THE POSITION OF ITS RIGID PERFECT INVOLUTE COUNTERPART. NOTE ESPECIALLY THE 
SIGN CONVENTION FOR POSITIVE p 1=1,2, 3. ABOVE CONFIGURATION IS THE MIRROR 
IMAGE OF FIG. 6. 
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USE OF THE GENERALIZED TRANSMISSION ERROR IN THE 
EQUATIONS OF MOTION OF GEAR SYSTEMS 

In order to write the equations of motion of a gear system, one needs 
mathematical expressions for all forces and moments acting on each gear and each 
pinion in the system. Such forces and moments include those arising from bearing and 
shafting displacements and those transmitted by the teeth of the mating pinion or 
gear. Mathematical expressions for the bearing and shafting forces and moments can 
be written down in a relatively straightforward manner. In the present section, a 
methodology is developed for incorporating the forces and moments transmitted by the 
teeth of pairs of meshing gears into the equations of motion of gear systems. The 
effects of the nonuniform transmission of motion between meshing gears arising from 
shafting translational and rotational displacements, tooth and gear body elastic 
deformations, and deviations of the tooth running surfaces from perfect involute 
surfaces are included in this methodology. 

For pairs of rigid perfect parallel-axis (refs. 7-9) and bevel involute gears, 
all tooth contact takes place in a fictitious plane called the plane of contact or 
base plane. If the forces of friction between the teeth are neglected, then the 
forces transmitted by the teeth lie normal to the tooth surfaces, and, therefore, in 
the plane of contact (ref. 7). It follows that the resultant total force transmitted 
by all teeth simultaneously in contact can be described by three scalar components - 
namely, two orthogonal force components that lie in the plane of contact and the 
moment of the resultant total transmitted force taken about a specified point in the 
plane of contact. These three components are the components of the generalized mesh 
force vector (Fig. 4). 

In studying the vibratory excitation arising from a pair of meshing gears, we 
must be concerned with the unsteady behavior of the forces transmitted by the 
teeth. These unsteady forces are dependent on the tooth and gear body elasticity and 
the dynamic parameters of the gear system (inertias, stiffnesses, etc.) and therefore 
must be regarded as unknowns. In this section, a relationship is derived that 
relates the three-component generalized mesh force vector that describes these 
unsteady forces, the three-component generalized transmission error vector (Fig. 5) 
that describes the nonuniform transmission of motion, a three-component vector that 
describes the geometric deviations of the tooth running surfaces from perfect 
involute surfaces, and an appropriately defined mesh stiffness matrix. Combining 
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this relationship with the expression (2.12) for the generalized transmission error 
vector yields a matrix equation which is interpreted as an equation of constraint 
that relates the unknown vibratory displacements of the two meshing gears, the 
unknown mesh forces, and the geometric deviations of the tooth running surfaces from 
perfect involute surfaces which are assumed known. The relationship between the 
generalized transmission error vector, the generalized mesh force vector, and the 
geometric deviations of the tooth running surfaces from perfect involute surfaces is 
worked out in detail for spiral bevel gears of nominal (spherical) involute design 
and completely arbitrary spiral angle (19, p. 241). Earlier analyses along the same 
lines utilizing the one-component transmission error (refs. 3-12, 20) have been 
carried out for spur gears in refs. 6 and 20 and for helical gears in refs. 7 and 9. 


Equations for Generalized Transmission Error in Terms of Mesh 
Loading and Geometric Tooth Surface Deviations 


We now derive a set of equations for the generalized transmission error of 
spiral bevel gears that involves the mesh loading components illustrated in Fig. 4 
and the geometric deviations of the tooth running surfaces from perfect spherical 
involute surfaces. This material is an extension to the three-component generalized 
transmission error for spiral bevel gears of a much simpler derivation (ref. 7, pp. 
1412, 1413) of the one-component transmission error for helical gears. 


Force Deformation Relations. In the previous section, it was pointed out with 

the aid of Fig. 3 and equation (2.6) that specification of the pair of coordinate 

values x,y determines a unique point of contact on the running surfaces of two mating 

teeth designated by index j. Let u^*^(x,y) denote the elastic deformation of the 

running surface of tooth j of gear (•) under such a contact point, (•) ■ 1 or 2. The 

direction of u^*^(x,y) is taken normal to the tooth surface and, therefore, in the 

plane of contact normal to the line of contact, as can be envisaged with the aid of 

Fig. 3; u^*^(x,y) is defined to be positive when it is "equivalent" to removal of 

material from the running surface of the tooth. Thus, u'*'(x,y) always is non- 

nj 

negative. Let u n ^(x,y) denote the sum of the elastic deformations of two mating teeth 
in contact - i.e. , 


, , 4 ( 1 ). . ( 2 ). . 

u (x,y) = u . (x,y) + u (x,y) 
nj nj nj 


(3.1) 
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and let K^’^(x,y) denote the local stiffness of tooth pair j per unit length of line 
of contact as defined in Appendix B, where subscript n again denotes that these 
quantities refer to deformations normal to the common point of tangency of the two 
mating teeth which is located by the coordinate values x,y. The superscript (•) on 
K^*\x,y) will be used to distinguish the different definitions of local tooth-pair 
stiffness that are established in Appendix B and used later in this section. Then 


dW n j(x,y) - K^(x,y) u nJ (x,y) dl 


(3.2) 


is the normal force increment transmitted by the differential length dl of the line of 
contact of tooth pair j as illustrated in Fig. 9. From Fig. 9, we can see that 


dl 


sec >l> b (y) dy ; 


(3.3) 


hence , 


dW n j( x ,y) 


(•) 


sec * b (y) K nj 


(x,y) u n ^(x,y) dy 


(3.4) 


The angle ij> b (y) is the base plane spiral angle (ref. 19, p. 241) illustrated in Fig. 9 
which can vary in a rather general, but specified, way as a function of y * r Q - r. 


Let dW n j(x,y) be resolved into tangential and radial components dW t ^(x,y) and 
dW^j(x,y) that are perpendicular and parallel, respectively, to the zone of contact 
centerline as shown in Fig. 9. From the figure, it is readily seen that the differen- 
tial tangential and radial forces are related to the differential normal force by 


dw tj - sin{J- [* b (y) - Xj(x,y)]} dW nj 
- cos[i|> b (y) - Xj(x,y)] dW^ 


(3.5) 


and 


dW rj = cos{* - [+ b (y) - Xj(x,y)]} dW nj 
■ sin[i|> b (y) - Xj(x,y)] dW , 


(3.6) 
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■V 



FIG. 9. RESOLUTION OF DIFFERENTIAL NORMAL FORCE INCREMENT dW .(x,y) INTO COMPONENTS 

nj 

dW (x,y) and dW^(x,y) PERPENDICULAR AND PARALLEL TO ZONE OF CONTACT 
CENTERLINE. FOR LOCATION OF ZONE OF CONTACT IN PLANE OF CONTACT, SEE FIG. 3. 
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where A^(x,y) is the angle in the base plane between the zone of contact centerline 
and the radial line passing through the differential force element with vertex at the 
base cone apex as illustrated in Fig. 9. The contributions W t j(x) and W r j(x) from the 
force transmitted by tooth pair j to the tangential and radial components W t ( x ) and 
W r (x) of the total force transmitted by the mesh are obtained by (Stieltjes) 
integrations of dW t j and dW r j , respectively, over the total line of contact of tooth 
pair j: 

y B 

w t j( x > = / dW tJ (x,y) (3.7) 


y B 

W ri (x) = / dw ri ( x >y) » (3.8) 

y A 

where y^ and yg denote the endpoints in the coordinate y of the line of contact of 
tooth pair j. 


Let h denote the distance measured perpendicular to dW n j( x ,y) between the zone of 
contact center 0^ and the line of action of the differential normal force dW n j(x,y) as 
illustrated in Fig. 9. Then, using the sign convention for moments illustrated in 



pair j to the total moment t(x) about the zone of contact center 0 is 


y B 

Tj(x) = - / h(x,y) dW n j(x,y) . (3.9) 

y A 


The tangential and radial components W t ( x ) and W r (x) of the total force transmitted by 


the mesh and the 
by the sums of W 

total moment x(x) about 
t j(x), W r j(x), and Tj(x) 

the point _0_ illustrated 
over all tooth pairs in 

in Fig. 4 
contact : 

are given 

W fc (x) = 

I 

j 

V 10 



(3.10) 

W f (x) = 

l 

3 

V x) 



(3.11) 

t(x) = 

l 

j 

TjCx) . 



(3.12) 
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Inserting equations (4) and (5) into equation (7) gives 


/>■) 

W (x) = / sec i|) fe (y) K nj ( x »y) cos[* b (y)-Aj(x,y)] u nJ (x,y)dy , (3.13) 

y A 

and inserting equations (4) and (6) into equation (8) gives 
78 ( A 

W^(x) “ / seo|/ b (y) (x,y) sin|> b (y)-Aj (x,y)] u n ^(x,y)dy , (3.14) 

y A 

where K^^(x,y) and K^^(x,y) are defined in Appendix B. In Appendix C, it is shown 
nj nj 

that 


h(x,y) = r o cos[i|> b (y)-Aj(x,y)] - (r o -y)cosi|» b (y) . (3.15) 

Finally, inserting equations (4) and (15) into equation (9) yields 


Tj(x) = -r Q / seo|i b (y) K^(x,y) cos[i|> b (y)-A.j (x,y)] u n ^(x,y)dy 


y B ( > 

+ J (r -y) K^ ; (x,y) u n (x,y)dy , 
y A 


(3.16) 


where K^^(x,y) also is defined in Appendix B. 
n J 

Force-Transmission Error Relations. Let q^*^(x,y), (0=1 or 2, denote the 

geometric deviation of the unloaded running surface of tooth j of gear (• ) from the 

surface of its perfect spherical involute counterpart at the point of contact on the 

surface determined by the pair of coordinate values x,y. The deviation ri^*^(x,y) is 

measured in a direction normal to the tooth surface and, therefore, in the plane of 

contact normal to the line of contact; q^*^(x,y) is taken to be positive when it is 

n J 

equivalent to removal of material from the surface of its perfect involute counter- 
part. Let r)^(x,y) denote the sum of the geometric deviations of two mating teeth in 
contact - i.e., 
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(3.17) 


T 1 .(x,y) = r/^(x,y) + i/^(x,y) , 
nj nj nj 

which is the geometric deviation counterpart to the corresponding relation ( 1 ) that 

applies to elastic deformations. Let £ . (x,y) denote the local transmission error of 

nj 

tooth pair j defined in the same direction and with the same sign convention as 
v^jCx.y) and n n j(x,y). The local transmission error is defined as the sum of the 
elastic deformation u n j(x,y) and the geometric deviation h n j(x,y) of the tooth pair at 
the point of contact determined by x,y: 

£ .(x,y) - u (x,y) + n (x,y) . (3.18) 

n J nj nj 

In Fig. 10, the local transmission error t .(x,y) of tooth pair j is shown 

nj 

resolved into local tangential and radial components £ t - ap and £^ - by that are 
perpendicular and parallel to the zone of contact centerline, respectively. The sign 
convention of the transmission error components shown in Fig. 10 coincides with the 
displacement of the gear on the left-hand side of the zone of contact shown in Fig. 5. 
The local tangential component £ - ay contains contributions from the global 

tangential and rotational components 5 and p that are illustrated in Fig. 5; the 
local radial component £^ - bp contains contributions from the global radial and 
rotational components £^ and p also shown in Fig. 5. From Fig. 11, one can readily 
show that the distances a and b (in Fig. 10) of the contact point from the center 0 _ of 
the zone of contact are 

a * r Q [ l-cosXj(x,y)] + y cosXj(x,y) (3.19) 

b ■ (r Q -y) sinXj(x,y) . (3.20) 

Thus, the local tangential and radial components of the transmission error are 

C fc (x) - ap(x) = £ fc (x) - {r o [l-cosXj(x,y)] + y cosXj(x,y)} p(x) (3.21) 

and 

C (x) - bp(x) = £ (x) - [(r -y) sinX.(x,y)] p(x) . (3.22) 

r r L o j J 

Recognizing the geometric similarity between the force resolution in Fig. 9 and the 
local transmission error resolution in Figs. 10 and 11, it follows from equations (5) 
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LINE OF CONTACT 
OF TOOTH PAIR j 


\ 





FIG. 10. RESOLUTION OF LOCAL TRANSMISSION ERROR C nj (x,y) INTO LOCAL TANGENTIAL AND 
RADIAL COMPONENTS ^-ay AND ? -by PERPENDICULAR AND PARALLEL TO ZONE OF 
CONTACT CENTERLINE, RESPECTIVELY. DISPLACEMENTS SHOWN COINCIDE WITH MOTION 
OF GEAR ON THE LEFT-HAND SIDE OF THE ZONE OF CONTACT SHOWN IN FIG. 5. 
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and (6) that the tangential and radial components of the local transmission error are 
related to the normal component by 


5 t (x) - ap(x) = cos[ij) b (y) - Xj(x,y)] 

and 

C r (x) - bp(x) = sin[r|) b (y) - Xj(x,y)] 


nj <*,y> 

(3.23) 

nj <x - y) • 

(3.24) 


It is possible, now, to express the (unknown) elastic deformation component 

u n j(x,y) in equations (13), (14), and (16) in terms of the geometric deviations 

q ,(x,y) of the tooth running surfaces and the three generalized transmission error 
*Vj 

components p(x), ? t (x), and e r (x). Solving equation (18) for u n j( x ,y) gives 

u nj (x,y) = C nJ (x.7) “ n nj (x,y) . (3.25) 

Multiplying equation (25) by cos[t(> b (y) - Xj(x,y)] and combining the resulting 
expression with equations (21) and (23) to eliminate ? n j(x,y) gi ves 

cos[i|» b (y) - Xj(x,y)] u^x.y) = ? t (x) - { rj 1 - cosXj(x.y)] + y cosXj(x,y)} u(x) 

- cos[^ b (y) - Xj (x,y)] n n j(x,y) . (3.26) 

Multiplying equation (25) by sin[t|» b (y) - Xj(x,y)] and combining the resulting 
expression with equations (22) and (24) to eliminate 5 n j(x,y) gi ves 

sin[i}> b (y) - Xj(x,y)] u n ^(x,y) = e r (x) - [ (r Q -y) sinXj(x,y)] y(x) 

- sin[^ b (y) - Xj(x,y)] h n j(x,y) . (3.27) 

Dividing equation (26) by cos[i|> b (y) - Xj(x,y)] gives 

u (x,y) = secU (y) - X . (x,y)l ft (x) - {r [ 1-cosX (x,y)] + y cosX (x,y)} u(x)) 
nj d j t o j j 

- q .(x,y) . (3.28) 

n J 

Finally, substituting equations (26), (27), and (28) into equations (16), (13), and 
(14), as required, yields the desired expressions for tj(x), (x) , and W^^x): 
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y B 

T.(x) = -r / sec <P b (y) K^j ) (x,y)(c t (x) - { rj 1-cosX, (x,y)] 
y A 

+ y cosXj (x,y)} y(x) - cos[i|» b (y) - Xj(x,y)] b n j(x,y)) dy 
y B 

+ / (r Q -y) K^ ) (x,y)[secti(» b (y)-X (x,y)](c t (x) 

y A 

- {r Q [ 1-cosX j(x,y) ] + y cosXj(x,y)} y(x)) - n n j(x,y)] dy , (3.29) 


W t j(x) - / sec 4> b (y) K^\x,y)(^ t (x) - { t q [ 1-cosX j (x,y)] + y cosXj(x.y)} y(x) 


- cos[* b (y) - Xj(x,y)] n n1 (x,y)) dy , 


nj 


and 


(3.30) 


W rj (x) - / sea|» b (y)K^ ) (x,y){; r (x) - [ (r Q -y)sinXj (x,y)] u(x) 


-sin[i|> b (y)-Xj(x,y)] Ti n ^(x,y)} dy. 


(3.31) 


Let us define 


y B 

(x;t|)^) - / sec 4> fe (y) K^(x,y) dy 




(3.32a) 


«(•) 


y B 


<•), 


K v ' ' (x;xp fe — X^ ) - / sec[i|i b (y) - Xj(x,y)] K r j (x,y) dy 




(3.32b) 


(x;iJ» b> X ) * / sec i|> b (y) sinXj(x,y) K^(x,y) dy 

c-i S j y A 


(3.32c) 
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* (xj* ,x ) =/ sec i(), (y)[ 1-cosX . (x,y)] K^ \x,y) dy 

cr-l(l-c)j y A J nj 


~(.) A B ( \ 

K (x;i)» -X ,X ) =/ sec[i|i b (y)-X.(x,y)][l-cosX .(x,y)] K’'* ; (x,y) dy ( 

c— 1 (l-c)j J y A J J 


;<•> 


K vw ( x; ^ b - X j ) = / y sec[i|) b (y) - Xj(x,y)] K^Cx.y) dy 


( 


ycr-ij 


~(*) A B 

K (xJ'I'.X ) = / y sec i|v(y) cosX.(x,y) K^*.\x,y) dy 

yc-lcj 3 y A nj 


C 


K (x;i(» -X ,X ) - J y sec[^ b (y)-X .(x,y)] cosX,(x,y) IT* ; (x,y) dy (: 

yc -1 c j 3 3 y A 3 3 


K (,) 

ycr-lsj 


A B (.) 

< x ;Vy “ / y ® e «l» b (y) sinX j(x,y) r* ; (x,y) dy 


c 


K 


(•) , .... A ' B 


yc~ 1 (l-c)j 


(x5 *b" X j ’V “ / y sec ['*'b (y) “ X j (x ’ y) H 1 ~ cosX j( x * y )] K nj^ x « y ) dy 


0 


K (x;i|/ -X ,X ) ■ / y 2 sec[^ b (y)-X.(x,y)]cosX.(x,y) K^*\x,y) dy (' 

y 2 c — 1 c j ° 3 3 y A J J J 


“ / K^^Cx.y) ti . (x ,y) dy 


Kj 


nj "nj 


(• 


5.32d) 
l.32e) 
1. 32f ) 
l.32g) 
i.32h) 
.321) 

• 32 j ) 
.32k) 

.320 
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~(. ) A r B (.) 

n (x;<J> K ,i|)^-X J = j K nj ^ x#y ^ sec ^ h (y) cos[i|) h (y)-X 1 (x,y)] n n1 (x,y) dy 


Kc-lcj 


b b j 


r h ,j, "nj 


(3.32m) 


#w ( i ) . / v 

n (x;iJ) b> i|) b -X^) = J K n * ; (x,y) sea|) h (y) sin[^ l> (y)-X < (x,y)] n n .j(x,y) dy 


Kcflsj 


V'' >l n j 


(3. 32n) 


Vj (x) " ^ y K nj )(x ’ y) n nj (x ’ y) dy • 


(3.32o) 


The above notational convention is easy to remember once it is understood. Each of 
the above fifteen definitions is an integral over the full line of contact of tooth 
pair j. The integration is represented by the tilde over the symbols on the left-hand 
side - since the tilde has the approximate form of an integral sign (rotated 90°). 

The fact that the integration is taken over the line of contact of tooth pair j is 
represented by the last of the subscripts, j, on the symbol the tilde has been placed 
over. The first eleven of the definitions are regarded as weighted integrals of the 
stiffness K^\x,y) and the remaining four are regarded as weighted integrals of the 
deviation Ti n j(x,y), which explains the symbols that the tildes have been placed 
over. The subscripts following these symbols represent the weighting functions in the 
integrals, where c — 1 represents sec = cos — 1 , s represents sin, c represents cos, 1-c 
represents 1-cos, and K represents K^*^(x,y). Each of the fifteen integrals is a 
function of x = r Q 0, which accounts for the symbol x in the arguments on the left-hand 
sides. The remaining symbols after the semicolons in the arguments on the left-hand 
sides represent the arguments of the subscripted trigonometric functions, in the same 
order as the subscripts, with the last argument on the left-hand side after the 
semicolon in each of the definitions representing the argument of the last subscripted 
function before the j. Each superscript (•) in the left-hand side represents the 
superscript (•) in K^*\x,y) in the right-hand side, (• ) = (t), 

(r), or (y), as appropriate. 

By utilizing the definitions (32a) - (32o), equations (29), (30), and (31) can be 
expressed as 
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t,(x) = -r r (x) K (t) <xj* ) + r 2 y(x) K (t) (x;i,L) 

j o C c~lj b ° c-l(l-c)j b J 


+r u (x) (x;i|> ,X.) + r (x;i|» ,\|> -X .) 

0 yc-lcj b j ° Kc-lcj b b j 


+r r (x) (x;i|i -X.) - ? (x) (x;tp -X ) 

° c c -ij b j t yc _ ld b j 


-r 2 y (x) K (u) (x;i().-X.,X.) + r u(x) K kv ' - -X . ,X , ) 

° b j j o l D J J 


(u) 


c 1 ( 1— c) j 


yc l(l-c)j 


-r y (x) (x;i|>, -X . ,X.) + y(x) (x ;\p, -X . ,X . ) 

0 yc-‘c] b 3 3 y 2 c — 1 c j b 3 3 


-r )(x) + n^(x) 
o n Kj ^ ' n yKj ' ’ 


(3.33) 


W ,(x) = c (x) (x;i|» ) - r y (x) (x;ij/ ,X.) 

tj * c-lj b ° c-l(l-c)j b i 


-y(x) (x;i|i.,X.) - ri^ (x;ij> -X.) , (3.34) 

yc^cj b j Kc-lcj b b 3 


and 


W .(x) * C_(x) (x;ij) ) - r y (x) (x;ij» ,X.) 

rj r c-lj b ° c _ lsj b 3 

+ y(x) (x;\J» ,X.) - (x;tj> ,t|; -X.) . (3.35) 

yc-lsj b j Kc-lsj b b 3 

The total moment t(x) and forces W fc (x) and W r (x) transmitted by the mesh are obtained 
by summing equations (33), (34), and (35) over all tooth pairs j in contact, as 
indicated by equations (10), (11), and (12). If the sums over j of the quantities 
defined by equations (32a) - (32o) are denoted by the same symbols with the subscripts 
j that are not included within the parentheses omitted - i.e., 
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K^(x;i|» ) = l (x;4> ) 

b . i . b 


(3.36a) 


n-1 


j c-lj 


\x;* -X ) = l \x;i|) -X .) 
c-1 b J j c-lj b J 


(3.36b) 


~(* ) , v A v ~(* ) . . 
n _ (x) = l n (x) , 
yK j yKj 


(3.36o) 


then the total moment and forces transmitted by the mesh are expressed in a form 
identical to equations (33), (34), and (35) with the subscripts omitted as indicated 
in equations (10) - (12) and (36a) - (36o). The coefficients of p(x), ? fc (x), and 
C r (x) in the resulting equations can be factored out to yield the following set of 
equations for the generalized transmission error components: 


y(x)(K (w) (x;^ -X ,X ) + r [K (t) (x;* X)-K (y) (x;i|) -X ,X )] 

v2 p lp J J ycT’l C ^ \7r* — 1 n ^ J J 


y<- C 1 C 


ycr l c 


+ •- + r „t* <t . ) . - k( “.\ _ .x^)] n 


v * Th A A 

yc - 1 (l-c) ^ ^ 


c-l(l-c) b J c-1 (l-c) b J J 


?t (x){K (w) (x;^ b -X ) + r [K (t) (x;tf. ) - K (y ) (x;^ -X )] } 
yc~ 1 J c“ 1 c ~ 1 J 


= t (x) - hy^(x) + r o [n^ y \x) - (x;i|i b ,i|) b -Xj )] 


(3.37) 


-y(x)[K^ (x;i|) ,X ) + r (x;tJ» ,X . )] + £ . (x)K^(x;\|> ) 

yc - 1 C b j ° r-m-rl b J t c-l b 


C 1 (l-c) 


W fc (x) + (x;ij> b ,ij> b -x .) 

Kc _ 1 c J 


(3.38) 
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and 


-y(x)[ r o K^ (x;i() b ,X ) - (xj^.X .)] + 5 r (x)K^(x;i|i ) 

c —1 s J yc — ■ l s J c — 1 

= W(x)+n (r) (x;ip,,il>-X.). (3.39) 

1 Kc-ls b b 2 

Equations (37) and (38) are a pair of simultaneous linear algebraic equations for 
the rotational and tangential components, y(x) and £ t (x), of the generalized tran- 
smission error illustrated in Fig. 5. The coefficients of the unknowns defined in 
equations (32a) - (32k) are functions of the gear mesh geometry and the tooth-pair 
stiffnesses K^*\x,y), (•) = (t) or (y ) . The right-hand sides contain the forcing 
terms x(x) and W fc (x) illustrated in Fig. 4 and the stiffness weighted geometric 
deviations of the tooth running surfaces from perfect involute surfaces defined in 
equations (32i) - (32o). Equation (39), the third in the set, can be regarded as an 
equation for the radial component £^(x) of the generalized transmission error in terms 
of y(x) and the forcing and geometric deviation terms that are written on its right- 
hand side. 


Use of Generalized Transmission Error in Equations of Motion 

The set of equations (37) - (39) for the generalized transmission error 
components y(x), ? t (x), and £ r (x) can be written in matrix form as 

K(x) £(x) = W(x) + n K (x) , (3.40) 

where 5 (x) is the generalized transmission error vector defined by equation (2.10), 
W(x) is the generalized mesh force vector 

t(x) 

W fc (x) 

W r (x) 

an< * '*' S tbe vector that describes the stiffness weighted geometric deviations of 

the tooth surfaces from perfect involute surfaces: 



(3.41) 
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-n^\x) + r [r/ P \x) 
yK 0 K 


~(t) 


n (t ^ (x;i|> ,i|> -X . )] 
Kc~ 1 c b b J 


n (x) m < n -x ) 


K 


Kc-1, 


b b j 


^v(r) 

n (x;<j» ,<|> -x.) 
(. KC - 1 s b b J 


(3.42) 


The matrix K(x) is the 3x3 element coefficient matrix readily written down from the 
left-hand sides of equations (37) - (39). 

Let us define the inverse C(x) of the mesh stiffness matrix K(x) by 

C(x) = K-l(x) . (3.43) 


An expression for C(x) is derived in Appendix D. Premultiplying equation (40) by 
K“ l (x) gives 


C(x) = C(x)W(x) + C(x)n K (x) , (3.44) 

where the definition (43) has been used. Equating the right-hand sides of equations 
(2.12) and (44) yields 

T (1) q (1) (x) + T (2) q (2) (x) = C(x)W(x) + C(x)n K (x) 
or 

T (1) q (1) (x) + T (2) q (2) (x) - C(x)W(x) - C(x)n K (x) , 

which is the final result. 

Each gear and pinion in a simple or complex gearing system possesses 6 degrees- 
of-freedom when modeled as a rigid body. These 6 degrees-of-freedom can be taken to 
be the x£* ^ and P^*^» i = 1>2,3 illustrated in Figs. 6 or 8, as appropriate. The 
differential equations of motion of each such gear or pinion can be expressed in terms 
of these 6 coordinates and their first and second temporal derivatives, the elastic 
restoring forces and moments of the bearing supports which depend on these same 6 
coordinates, and the shafting forces and moments that also may depend on the compar- 
able coordinates that describe the position of a neighboring gear or pinion. T.n 


(3.45a) 

(3.45b) 
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addition to these elastic forces, the three components of the mesh force vector 
W, equation (41) and Fig. 4, also act on the gear. The temporal variations of the 
components of the mesh force vector W depend on the dynamic properties of the gearing 
system (masses, stiffnesses, etc.); therefore, they must be treated as dependent 
variables or unknowns whose time-dependent behavior also must be obtained from solu- 
tion of the equations of motion of the gearing system. 

Treatment of the three components of W as unknowns requires the addition of three 
equations to the equations of motion of the gearing system for each pair of meshing 
gears in the system. The three scalar components of the matrix equation (45b) are 
these three additional equations. If, for each pair of meshing gears in the system, 
we consider the components of W as three additional generalized coordinates (ref. 21, 
p. 12), then the scalar components of equation (45b) have the form of holonomic 
constraints (ref. 21, pp. 10-14; ref. 22, pp. 24-27). When a matrix equation of the 
form of equation (45b) is included for each pair of meshing gears in the system, it is 
evident that the vibratory excitation of the system arising from each gear pair is a 
displacement form of excitation provided by the three components of the vector 
C(x)tj^(x) together with the three components of the vector C(x)W(x) that provide a 
parametric excitation of the system through the periodic dependence on x of the matrix 
C(x). 


Denoting the nominal angular velocity of gear (•) by rad/sec, it follows 

from equation (2.8) that 

x - r£;> (e (,) - e^) * r£;> 0 ) (,) (t - t o ) , (3.46) 

where 9^*^ and t denote appropriate initial conditions. From equation (46) it 
follows that C(x) can be regarded as a function of time t. Hence, the equation of 
constraint (45b) has time-dependent coefficients which sometimes are referred to as 
rheonomic constraints (ref. 22, p. 32). When included as part of the equations of 
motion of a gearing system, it seems appropriate to refer to the matrix equation (45b) 
as the mesh constraint equation. However, from a physical perspective, equation (45b) 
is a geometric compatibility relation (ref. 23, pp. 99, 100) involving the gear axes 
vectorial displacements q^^(x) and q^^(x), the elastic deformations of the gear 
teeth and bodies characterized by the matrix product C(x)W(x), and the geometric 
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deviations of the tooth running surfaces from perfect involute surfaces characterized 
by the matrix product C(x)n^(x) - as is readily understood from equations (2.11) and 
(2.12), and equations (44) and (45) above. 
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EFFECTS OF BEARING OFFSET AND FLEXIBILITY ON THE 
MESH FORCE DISTRIBUTION 


In the preceding section, a method is outlined for determining the values of the 
three generalized force components W t , w r , and x. This method is valid for elastic 
teeth of modified spherical involute design if the true surface of tooth contact is 
reasonably well approximated by a plane surface passing through the common apex of 
the base and pitch cones of the meshing bevel gears. The method outlined in the text 
following equation (3.45b) includes contributions to W t> w r , and x from the inertial 
forces of the gear bodies arising from their transverse and axial vibratory motions , 
and contributions from the torsional vibratory motions of the gears which depend on 
the torsional dynamic properties of all gearing elements in the system. 

In the present section, equations are derived for the generalized mesh force 
components W t , W r , and x under the assumption that the gear body inertial forces and 
the inertial forces associated with the "effective" bearing mass are negligible in 
comparison with the bearing/bearing support elastic reaction forces. When the values 
of these generalized mesh force components are used in equations (3.41) and (3.44), 
the value of the generalized transmission error determined therefrom as a function of 
the rotational positions of the gears is the generalized static transmission error. 
However, contributions from the torsional vibratory motions of the gears are readily 
taken into account when using the equations derived herein, since the component of 
mesh torque coinciding with the gear axis is an "input" to the equations for W t , W r > 
and x . The temporal variations of this torque component can be determined from the 
torsional equations of motion of the gearing system. 

In the subsection below, the generalized mesh force vector is transformed into a 
coordinate system suitable for use in writing the equations of motion of the gear 
bodies. This transformation is useful in applying the methodology outlined following 
equation (3.45b) which takes into full account the inertial forces associated with 
transverse and axial gear body motions. 

Expression for Mesh Forces in Gear Coordinates 

The scalar components W^, W r > and x of the resultant total force W transmitted 
by the gear teeth can be represented by the generalized mesh force column matrix or 
vector defined by equation (3.41) and illustrated in Fig. 4, 
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w = < w 


(4.1) 


According to Newton's third law, this system of forces acts on each member of the 
pair of meshing spiral bevel gears. 

We require a representation in “gear coordinates" of the mesh force system W fc , 

W , and x acting on a generic spiral bevel gear (•). These gear coordinates are 
r \ 
provided by the generalized displacement column matrix or vector defined by 

equation (2.11), 


(• ) A . 

q -■< 


P3 


(•> ^ 

(•) 

(•) 

1 

(•) 
(•) 
(•) 




(4.2) 


where each p i = 1,2,3 represents the (small) rotational deviation about the 

i (.) 

base cone coordinate axis x^ of the position of the shaft centerline of gear (. ) 

from the position of its rigid perfect involute counterpart. Orthogonal coordinates 

(.) (.) (.) 

x^ * x 2 * and x 3 represent translations of the shaft centerline position of 


gear (•) from the position of its rigid perfect involute counterpart. This six 
degree-of-freedom system of coordinates is illustrated in Fig. 12 or its mirror 
image Fig. 13, whichever is appropriate. 

We can associate a generalized force column matrix or vector 


r (•) > 

T 1 

(•) 


q ( ‘ ) a H 


^ F 3 


(•) 

l 

(•) 

(•) 

<•) 


(4.3) 
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FIG. 12. GEAR COORDINATE SYSTEM CONSISTING OF ROTATIONS p i = 1,2,3 AND 

TRANSLATIONS x^’^, i - 1,2,3 REPRESENTING DEVIATIONS OF POSITION OF GEAR 

(•) FROM POSITION OF ITS RIGID PERFECT INVOLUTE COUNTERPART. CARTESIAN 

COORDINATES z ± , 1 = 1,2,3 ARE ORIENTED WITH RESPECT TO PLANE OF CONTACT 

WHICH IS z x , z 3 PLANE WITH z 3 AXIS PASSING THROUGH ZONE OF CONTACT CENTER 

_0. FORCE SYSTEM LOCATED AT 0_ SHOWS FORCES EXERTED ON GEAR (•) BY ITS 

MATING GEAR WITH MOMENT t OF RESULTANT TOTAL FORCE TAKEN ABOUT 0. FORCE 

SYSTEM LOCATED AT BASE CONE APEX SHOWS EQUIVALENT FORCES WITH MOMENT 

t> + r W OF RESULTANT TOTAL FORCE TAKEN ABOUT BASE CONE APEX. VECTOR 
0 t 

REPRESENTATION OF MOMENTS t and t + r W IS GOVERNED BY RIGHT-HAND RULE. 

o t 
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FIG. 13. MIRROR IMAGE OF FIG. 12. SEE CAPTION OF FIG. 12. HERE, VECTOR 


REPRESENTATION OF MOMENTS t and t + r W IS GOVERNED BY LEFT-HAND RULE 


with the generalized displacement vector q^* ^ given by equation (2). The general- 
ized force vector has three moment or torque components x^ \ i = 1,2,3 and 

(. ) i 

three force components , i = 1,2,3. These components are not the generalized 

forces normally associated with Lagrange's equations. 


Each i ■ 1,2,3 is the component of torque taken about the axis x^* \ 

1 (. ) i 

Positive values of the x^ correspond to positive values of the (small) rotations 

\ The directions of positive p^* \ = 1*2,3 are illustrated in Fig. 12 or 

Fig. 13, as appropriate. From Figs. 12 or 13, we may conclude that each torque 

(. ) 

component x^ is taken about the point occupied by the base cone apex. From Fig. 

12, one may observe that for each i = 1,2,3 the relationship between the sign of 

(. ) (. ) (. ) 

p^ , and hence x^ , and the sign of x^ is governed by the well-known right- 

hand rule; whereas, from Fig. 13, which is the mirror image of Fig. 12, the 

(.) (.) (.) 

relationship between the sign of p^ , and hence x^ , and the sign of x^ is 

governed by the left-hand rule. This difference in sign convention between Figs. 12 

and 13 is a necessary consequence of the mirror image property. 

(.) 

Each F , i = 1,2,3 in equation (3) is the component of force in the direction 

(.) i (.) 

of x^ . Positive values of the F^ correspond to positive values of the 
displacements x ^ 


For a generic spiral bevel gear (•), it is shown in Appendix E that the 
generalized force vector in gear coordinates is related to the generalized mesh 

force vector W by the matrix transformation 


.(•> 


Q' ' - -T 'w , (4.4) 

where ^ is the transpose of the matrix T^* ^ defined by equation (2.13) and Table 1, 


:(•) A 


0 




(•) 


“*Y 


(•) 


0 


"o^b 


-r sy 
o 'b 


-cv 


(•) 
(•) 

(•) 


(•) (•) 
-SY SV 

b 

(•) (•) 
-cy sv 

b 


± sv ( -> 


(•) <•) 

+ sy cv 

b 

_ (•) (•) 
+ cy cv 

b 


(4.5) 
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where c and s are abbreviations (ref. 18, p. 106) for cosine and sine operations, 

respectively, on the angles an< i These angles and r Q are illustrated in 

Figs. 12 and 13. Determination of the appropriate sign for the terms in the third 
-(•) 

column of T is governed by the rules for determination of the appropriate sign for 
the terms in the third row of T^* ^ that are described in Table 1. Further explana- 
tion of the gearing configurations described in Table 1 may be found following 
equation (2.13). 

The transformation provided by equations (4) and (5) yields the three torque 

components x.^*^ and three force components F.^'\ i « 1,2,3 associated with the 
1 / \ 1 / \ 

three rotations *' and three displacements x^* , i * 1,2,3 with positive torques 
and forces associated with positive rotations and displacements, respectively. These 
are the torque and force components applied to gear (») by its mating gear expressed 
in the gear coordinate system defined by equation (2) and Figs. 12 or 13, as 
appropriate. 


Expression for Generalized Transmission Error in Terms of Bearing Offset 

and Flexibility Matrices 

The expression given by equation (2.12), 


C (x) = T ( 1 ) q U) (x) 


+ T ( 2 Y 2 >(x) 


(4.6) 


describes the generalized transmission error $(x) of a meshing pair of spiral bevel 
gears in terms of the generalized displacement vectors q^^(x) and q^^(x) of the 
shaft centerline positions of gears (1) and (2), respectively, measured from the 
positions of their rigid perfect involute counterparts. The generalized transmission 
error is defined by equation (2.10) and Fig. 5; the matrix T^*\ (•) * (1) or (2), is 
defined by equation (2.13) and Table 1; and, the vectors q^*\ (•) = (1) or (2), are 
defined by equation (2). 

Geometric compatibility requires that the shaft centerline positions of gears 
( 1) and (2) coincide with their respective bearing centerline positions. The bearing 
centerline positions also are described using the generalized displacement vectors 
q^^(x) and q^^(x) that characterize the deviations of the actual centerline 
positions from the positions of the rigid perfect involute counterparts to the gears 
under consideration. These deviations have three vector contributions. The first of 
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(.) 

these contributions is the generalized displacement vector , (• ) * (1) or (2), 

that describes the deviation of the bearing centerline position when the gears are 
transmitting negligible load. This bearing "offset'' from the position occupied by 
the shaft centerline of the rigid perfect involute counterpart may be intentional or 


unintentional. The second contribution 


(•) 


describes the component attributable to 


the bearing and bearing support elastic deformations caused by the generalized force 

W transmitted by the gear mesh. Neither of the vector contributions or q ^ \ 

(.)(.) 0 " 

— 4 1 .. J t- - v ' __ _ > / _ c i 


by definition, includes a nonzero component q 


of the generalized 


/ \ 

displacement vector q v *' that describes the pure rotational deviation of the shaft of 
gear (• ) about its axis from the position of its rigid perfect involute counter- 
part. The vector describing this third contribution is {6 p ^ where 6.- is 

13 * 3 ^ 

Kronecker's delta (ref. 24). Hence, the third component of each vector 

q^ ^ and q^* ^ is, by definition, zero . The deviation vector for gear 
(• ) that describes the sum of these three contributions is 


(•), 

P- 1 * 


(4.7) 


(•) (*) ^ (O ,, 

I -0. * C if j 

Let ^ denote the 6x6 element bearing/bearing support flexibility matrix 


(•) 


[A 

contribution to the generalized displacement component q 


. ], where each flexibility influence coefficient A represents 

J m \ 

1 1 ^ ^ nn #■ n fVifl cranarol 4 vaA A A on 1 nn/imanf* nnmnnnnn 1 - « ' ' a nn-t h 


the 


from a unit torque or 


force component Q 


j 


(•) 


of the generalized force vector described by equations (3) and 


(4). Since the component q J" p ^ °f <1^* ^ is being handled separately by the 

/ \ 3 o . v 

term f‘l3»3 '1 in equation (7), the third row in the flexibility matrix *' is, by 

definition, filled with zeros; i.e.. 


(• ) A 
A = 0 
3j 


j = 1 , 2,.. .,6 • 


(4.8) 


Furthermore, it will be assumed that the reaction of the shaft of gear (•) to the 

(•) (•) 


pure torque component Q 3 
implies that 


is not coupled to the bearing supports, which 


(.) A 
A =0, 
i3 


1 , 2 ,... ,6 


(4.9) 


,<•> 

3 3 


since A„\' “ 0 by virtue of equation (8). Thus, the bearing/bearing support flexi- 


bility matrix [a|*^] is symmetric (ref. 23, pp. 4,5). The subscript • on 


(•) 
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Ls used as a reminder of the requirements (8) and (9). From the definitions of 
q^* ^ and A^ \ there follows 


(•) 




-v 


>~(< 

T 


(4.10a) 

(4.10b) 


where Q^* ' is defined by equation (3), and equation (4) has been used in going from 

equation (10a) to equation (10b). Equation (10b) gives the vector contribution of 

the deviation of the bearing centerline position attributable to the generalized mesh 

force vector W. From equations (9), (10a), and (3) it follows that q ^ 

contains no contribution from the component Q 3 “ T 3 of torque about the axis of 

gear (•), which may be determined from equations (4), (5), and (1) to be 

T~ ^ ^ = (x+r W )siny ^ \ 

3 v o t J b 

For (•) = (1) or (2), let us define the matrix 


(•) 4 _(• 


V 


)-(■ 

T 


(4.11) 


and the matrix 


# 


A (1) . (2) 

= * + # 


(4.12) 


Using equation (10b), we have 


T (l) S (1 > + - -[T (1 \ (1 ¥ 1) 4- T (2) A 0 (2) T (2) ]W 


= -* W , 


(4.13) 


where the definitions (11) and (12) have been used. Combining equations (6), (7), 
and (13) gives 


- «*V° ♦ I (2> ., 0 <2) - .»<*> ♦ T (,, ( Sl3 „ 3 (1) Wl ♦ , 


i3 


(4.14) 


where dependence on the rotational positions of the gears is denoted by the variable 
x. Equation (14) expresses the generalized transmission error ^(x) defined by 
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equation ( 2 . 10 ) in terms of the bearing centerline offset vectors <r_ and q n , 

^ (1) u 

the generalized mesh force vector W(x), and the gear shaft rotations p (x) and 

( 2 ) 3 

P 3 (x) about their axes illustrated in Figs. 12 or 13, as appropriate. 

The matrix * in equation (14) characterizes the composite bearing/bearing 

support flexibility of gears (1) and (2). Each of the two matrices and 

( 2 ) 

* that contributes to * in equation ( 12 ) is the negative of the bearing/bearing 
support flexibility matrix that yields the contribution to the generalized 
transmission error vector 5 from the bearing/bearing support elastic displacements 
caused by the generalized mesh force vector W. It follows by forming the transpose 
of ^ that each of these matrices is symmetric; i.e.. 


* = ( T Ao T ) “ T ) (T ) 


(T<-y 


(4.15) 


hence, it follows from the symmetry of ^ and the definitions ( 11 ) and ( 12 ) that 
^ and • are symmetric. Each of the matrices and * has three rows and 

three columns. 


Equations for Generalized Mesh Force Components 

Consider the last two terms in the right-hand side of equation (14). Utilizing 
equation (2.13) to evaluate these two terms, we have 


_,( 1 ) f . (1)/ u A 2), (2). 

{ 6 i3 P 3 (x)} + T {6 i3 p 3 (x)} 

( 1 ) ( 1 ). . ( 2 ) ( 2 ). . 

-sy b p 3 (x) -sy b p 3 (x) 


(1) (1), v (2) (2) 

- r sv P 3 00 -r *t P 3 00, 


(4.16a) 


( 1 ) 


(1) . (2), 


( 2 ). 


-[p 3 Vi/ (x)sinY b V ^ + p 3 v ^(x)sinY b VW ] 



(4.16b) 


54 



where, in going from equation (16b) to equation (16c), we have defined 


p (x) = p ^(x)siny ^ + p ^ 2 \x)sinY ^ . (4.17) 

3b 3 b 3 b 

Figure 14 shows the torque component x + r W at the base cone apex of gear 

/ v® ^ / \ 

(• ) resolved into components (x + r W )siny,^* ^ and (x + r W )cosYv^*' parallel to 

rv otb OtD 

the x^ ' and -x^ ' axes, respectively. In the absence of gear rotational inertial 

effects, it follows that the torque applied to gear (•) by its shaft, which coincides 
with the x q ^ ^ axis, must have a magnitude of (x + r W )siny, ^ ^ positive in the 
direction of — pg * . Examination of Figs. 12 and 13 shows this conclusion to be 
valid for all gear configurations listed in Table 1. Let us denote by ^ the value 
of this torque applied to the shaft of gear (• ); i.e., 

M^* ^ ■ (x + r W )siny ^ ^ , (4.18) 

o t b 

which we may regard as being known. Then, from equation (18), we have 

t * ^cscy ^ ^ - r W , (4.19) 

b o t 

from which it follows that the generalized mesh force vector, equation (1), can be 
expressed as 


>(•) (•) / \ 
M v ^cscy - r W (x) 
1 b o t 


W(x) = 


W^(x) 


W (x) 
r 


(4.20) 


where the dependence of W(x), W t (x), and W r (x) on the rotational positions of the 
gears has been explicitly denoted by their dependence on x. 

Equation 14 expresses the generalized transmission error c(x) in terms of bear 
ing offsets, bearing/bearing support flexibilities, and shaft rotational displace- 
ments. In contrast, equation (3.40), 




FIG. 14. 


GENERALIZED MESH FORCE COMPONENTS W t , W r , AND x SHOWN IN FIG. 4 EXERTED ON 
GEAR (.) BY ITS MATING GEAR. GEAR (• ) IS GEAR (1) IN FIG. 3. RIGHT-HAND 
RULE APPLIES TO VECTOR REPRESENTATION OF MOMENTS. MOMENT OF ft ABOUT ZONE 


OF CONTACT CENTER 0 IS t. MOMENT OF W ABOUT BASE CONE APEX IS x + r W 

WHICH HAS COMPONENTS (x + r W )sinY v ( * J AND (x + r W k )cosY v ( * ) IN THE 
(.) (.) o t b o t b 

x 3 AND -x DIRECTIONS, RESPECTIVELY. 
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(4.21) 


K(x)c(x) = W(x) + H R (x) , 


is an equation for the generalized transmission error in terms of the stiffness 
weighted geometric deviations n^x) of the tooth surfaces from perfect spherical 
involute surfaces defined by equation (3.42), and the tooth pair/gear body stiffness 
matrix K(x) defined by the coefficients of the left-hand sides of equations (3.37) 
through (3.39). By combining equations (14) and (21) to eliminate e(x), we can 
obtain a matrix algebraic equation for the generalized mesh force vector W(x). 


Substituting equation (16c) into equation (14) and the resulting expression into 
equation (21), we obtain 


K(x) 


C 3b <x) 

T U> q„ a> + T < 2 > q„ <2) - MU) - Jr p U) 

o Jb 
0 


= W(x) + h R (x) , 


(4.22) 


or, by rearranging terms and introducing the identity matrix I, 


c 3b (x) 

[K(x). + l]W(x) + K(x) ( r 0 p 3b ^ x V “ KUJfT^q,^ + T^q,^} - q^Cx) 

0 


(4.23) 


Introducing equation (20) into equation (23) and again rearranging terms, we have 


'~T W (x) s 
o t 


P 3b <X> 


[K(x)$ + I] <{ w^(x) \ + K(x) \ r 0 P 3b ( x ^ 


W (x) 
r 


K(x){T U) q 0 (1) + T (2) q 0 (2) } - h r (x) - [K(x)*+l]- 


„(•) (•)> 

M cscy, 

b 


0 

0 


(4.24) 
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Equation (24) is the desired matrix equation for the components W t (x) and W r (x) of 
the generalized mesh force vector W(x) defined by equation (1). Once W t (x) and W r (x) 
are obtained, the component t(x) of W(x) may be determined from equation (19) using 
the shaft input torque which is assumed known. 

The matrix equation (24) represents three simultaneous linear algebraic equa- 
tions for the three unknowns W t (x), W r (x), and p^^x). The (small) angular dis- 
placement vector component p^^(x) defined by equation (17) is the total contribution 
from both meshing gears to the rotational vector component normal to the base plane 
(plane of contact) arising from the angular deviation components and P 3 ^^ of 

the shafts of gears (1) and (2), respectively, from the positions of their rigid 
perfect involute counterparts. The value of p^(x) is of interest in the present 
investigation only in the supplementary role it plays in equation (24) as a 
contribution to our ability to compute W t (x) and W r (x). The parameter r Q is defined 
by equation (2.5) and the neighboring text. 

The matrix K(x) is the composite tooth pair/gear body stiffness matrix of the 
two meshing gears defined by the coefficients of the left-hand sides of the set of 
simultaneous equations (3.37) through (3.39), as one may see by comparing this set of 
equations with equation (3.40). The bearing/bearing support flexibility matrix 
defined by equations (11) and (12) occurs in equation (24) only as the second term in 
the matrix product K(x)# which collectively characterizes the elastic properties of 
the meshing gears (including the teeth) and the bearings and bearing supports. The 
matrix K(x), the matrix product K(x)*, and the identity matrix I all are of dimension 
3x3. 


The right-hand side of the matrix equation (24) contains contributions from 

three sources: (a) the bearing offset generalized displacement vectors q_^^(x) and 
( 2 ) 0 
q fl (x) from gears (1) and (2), respectively, (b) the vector q^(x) defined by 

equation (3.42) which characterizes the stiffness weighted geometric deviations of 

the tooth running surfaces from their perfect spherical involute counterparts, and 

(c) the torque * ^cscy^* \ where may be taken as the driving torque of either 

gear (•) of the meshing pair when the rotational inertial effects of gear (•) are 

negligible, and is the angle associated with the same gear (•) illustrated in 

D (1) 

either Fig. 12 or Fig. 13, as appropriate. The bearing offset vectors q (x) and 

( 2 ) ® 
q_ (x) are represented in the generalized displacement coordinates defined by 

equation (2) and Fig. 12 or Fig. 13, as appropriate. The matrix T'*' is defined by 
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SW f • } 

equation (2.13) and Table 1. It is the transpose of the matrix T defined by 
equation (5). Equation (24) readily can be solved using Cramer's rule (ref. 24). 

The location of the force vector W representing the resultant total force 
transmitted by the mesh is illustrated in Fig. 4. Dimension d shown in the figure 
measures the distance between the nominal zone of contact center _0_ located by r Q in 
Figs. 12 or 13 and the path of W. The value of d is readily computed from the three 
components x, W t , and W r of the generalized mesh force W described by equation (1). 
Since x = | W j d , it follows with the aid of equation (19) that 


d 



T 




(O (•) „ , . 

M cscy - r W (x) 
b o t 

1/2 

[W 2 (x) + W 2 (x)] 


(4.25) 


where W fc (x) and W r (x) are solutions to equation (24). Since these solutions are 
functions of the rotational positions of the gears, as indicated by their dependence 
on x, it follows that the distance d also is a function of the rotational positions 
of the gears. 


Case of Nonnegligible Gear Body Inertias. In deriving equation (24), the 
inertial forces arising from transverse and axial vibrations of the gear bodies were 
assumed to be negligible. This assumption was employed when the forces applied to 
the bearings by the gear bodies (or shafts) were set equal to the appropriate mesh 
forces; that is, when the generalized force vector Q^* ^ describing the forces applied 
to the bearings by the gear bodies (or shafts) was set equal to the right-hand side 
of equation (4). That step occurred in passing from equation (10a) to equation 
(10b). An assumption of negligible gear train inertia also was employed when the 
shaft input torque equation (18), was assumed to be known and independent of 

the rotational positions of the gears and of time as in equation (20). In addition, 
characterizing the relationship between the generalized force vector Q^* ^ applied to 
the bearings and the generalized bearing displacement vector q ^ by the bearing/ 
bearing support flexibility matrix = [A ] , as in equation (10a), ignores 

the inertia of the bearings and their supporting structure. When inertial effects 
are not negligible, a further complication arises from the fact that the mesh force 
components x, W t , and W r will be functions not only of the rotational positions of 
the gears characterized by the variable x, but also of time and the rotational speed 
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of the gear train which controls the temporal behavior of the gear train vibratory 
excitation. 


A general methodology for computing the components of the mesh force vector W in 
cases where inertial effects are not negligible is outlined in the text following 
equation (3.45b). Equations (4) and (5), which describe the components, equation 
(3), of the generalized force vector Q^* ^ in terms of the mesh force vector W, are 
useful for describing the mesh excitation torques and forces on the gear bodies in 
these cases. Furthermore, in these cases one would describe the bearing/bearing 
support elastic deformations by 



(4.26) 


where the meanings of 



(•) 

% ’ 



are the same as those in equation (7). 


Equation (24) remains valid if the input torque M^* ^ defined by equation (18) is 
allowed to vary with time. In such applications, M^* ^ is to be interpreted as the 
component of the mesh torque coinciding with the axis of gear (•). The temporal 
variations of M^* ^ can be determined by solving the torsional equations of motion of 
the gear system. If the temporal variations of M^* \correctly determined in this 
manner, are used in equation (24), the time dependent values of W t and W obtained 
therefrom, and the values of j obtained from M '* ' and W fc using equation (19), will 
include the inertial effects associated with torsional vibratory motions of the gear 
system. It is likely that these inertial effects are more important than those 
arising from transverse and axial vibrations of the gear bodies, or from bearing/ 
bearing support motions. 
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AN EXTREMUM PRINCIPLE FOR COMPUTATION OF THE ZONE OF TOOTH CONTACT 
AND GENERALIZED TRANSMISSION ERROR 

A set of linear algebraic equations, equations (4.24) and (4.19), has been 
derived for computing the three mesh force components under "static" loading 
conditions where inertial effects are negligible. The resulting mesh force 
components are dependent on the bearing/bearing support flexibilities, tooth/gear 
body stiffnesses, and other factors. When these mesh force components are used to 
evaluate the three transmission error components , they yield the generalized static 
transmission error components. 

To carry out the above-described computations, the positions of the endpoints of 
the lines of tooth contact must be known for all rotational positions of the gears. 
However, for tooth running surfaces that have been substantially modified from the 
"ideal” spherical involute toothform (ref. 13), these endpoint positions generally 
will have to be computed. 

In the present section, a method is developed for computing the positions of the 
endpoints of the lines of tooth contact and the three components of the generalized 
transmission error. For a given set of forces transmitted by the mesh, it is shown 
first that each of the three components of the generalized transmission error is 
stationary (ref. 22) with respect to small independent variations in the endpoint 
positions of all lines of tooth contact about their true values. For given forces 
transmitted by the mesh, this property implies that the value of each of the 
generalized transmission error components is insensitive to small errors in the 
endpoint positions of the lines of tooth contact. 

It then is shown that the tangential component of the generalized transmission 
error takes on a minimum value when the endpoint positions are given their true 
values. This property is the basis of a computational procedure utilizing the method 
of steepest descent (ref. 25) that is described for computing the true values of the 
endpoint positions and the generalized transmission error components. Alternative 
methods for computing the contact regions of the teeth of parallel-axis gears have 
been developed by Seager (ref. 26) and Conry and Seireg (refs. 27, 28). 

The computational procedure described herein is based on the assumption that 
along every path of tooth contact on a tooth running surface, the tooth surface 
modification from a perfect spherical involute surface is a monotonic strictly 
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increasing function (ref. 29) in both directions along this path away from some 
interior unmodified region which may be of negligible size. The modification is 
considered positive when it is equivalent to removal of material from a perfect 
spherical involute surface. However, one tooth of a meshing pair of teeth is 
permitted to be completely unmodified. If both teeth are unmodified, the zone of 
contact is the ideal geometric zone of tooth contact shown in Fig. 3. It will be 
assumed that the tooth running surface modifications from equispaced perfect spheri- 
cal involute surfaces are sufficiently small so that the true surface of tooth 
contact is well approximated by a plane surface passing through the common location 
of the apexes of the pitch cones of the two meshing gears as illustrated in Fig. 3. 

Stationary Property of Generalized Transmission Error 

As the gears in a meshing spiral bevel pair rotate, the lines of contact of 
meshing pairs of teeth pass through the nominal zone of contact, as one can see from 
Fig. 3. The variable x, defined by equation (2.6) and Fig. 3, describes the 
rotational positions of the gears and the locations of the lines of tooth contact in 
the plane of contact. Our aim is to determine the true locations of the endpoints of 
these lines of contact for a generic fixed value of x. When these endpoint locations 
are determined for all tooth pairs j in contact fop all values of x within the range 
“(A/2) < x < (a/ 2) where A is the tooth spacing interval in the variable x, the locus 
of points traced out by the endpoints of the lines of contact of all tooth pairs j as 
x is varied determines the zone of tooth contact in the plane of contact. This zone 
of contact is illustrated in Fig. 15. For a fixed value of x, the positions of the 
endpoints for tooth pair j are described by the values y^j and ygj where y is the 
radial coordinate 


y = r - r 
o 


(5.1) 


illustrated in Fig. 3. 


For a given generic value of x, it is shown below that each of the components 
? t (x), y(x), and C r ( x ) of the generalized transmission error is stationary with 
respect to small independent variations in the values of the endpoints y^. and ygj of 
all lines j of tooth contact of all pairs of teeth in contact, where these variations 
are taken about the true values of y^j or yg j . In carrying out these variations in 
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OF TOOTH PAIR] 


FIG. 15. ZONE OF TOOTH CONTACT IN PLANE OF CONTACT TRACED OUT BY ENDPOINT POSITIONS 
y A j AND y B j OF LINES OF TOOTH CONTACT. 


63 



the endpoint positions y^j or Yg j » all other parameters and variables, including the 
generalized force components W fc , t, and W r transmitted by the mesh, are held fixed. 

Furthermore, in carrying out these variations, negative values of the elastic 
deformations u^Cxjy) normal to the tooth surfaces are permitted, where a positive 
value of u n j(x,y) is associated with compressive forces between the teeth. The 
mathematical artifact of allowing negative values of u^jCx.y) to occur along the 
lines of contact outside the actual contact range determined by the true endpoint 
values of y^j and ygj will be seen later in this section to guarantee that the 
transmission error component £ fc (x) takes on a strong local minimum (ref. 30) with 
respect to both positive and negative variations in the endpoint positions about 
their true values. The procedure developed to compute the true values of the 
endpoint positions and the transmission error component C fc (x) also allows negative 
values of UjjjCx.y) beyond the true range of tooth contact. 

Stationary Property of c^Cx). Equations (3.37) and (3.38) are a pair of 
simultaneous algebraic equations for the generalized transmission error components 
y(x) and ? fc (x). These equations are of the form 

a U + at = r (5.2) 

I L t a 

b i y + V* = r v » (5.3) 

1 Z t D 

where the right-hand sides are of the form 


x + n. 


(5.4) 


w t + ,1 b ’ 


(5.5) 


where all quantities in equations (2) through (5) are functions of the variable x 

which describes the rotational position of the gears. The terms t and W t in 

equations (4) and (5), respectively, are two of the components of the generalized 

force transmitted by the mesh. Expressions for the terms a , a , n ; and b , b , 

12 a 1 2 

and Tig are readily obtained by comparing equations (2) through (5) with equations 
(3.37) and (3.38). Each of these latter six terms is an explicit function of the 
endpoint positions y^j and y B j of all of the tooth pairs j in contact at the 
rotational position of the gears designated by the value of x, as one can see from 
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equations (3.37) and (3.38), (3.36a) through (3.36o), and (3.32a) through (3.32o), 
where dependence of the endpoint positions y^j and y B j on tooth pair index j has been 
suppressed. Thus, the generalized transmission error components p and 
5 , determined by the solution of equations (2) and (3), also must be regarded as 
functions of the endpoint positions of all tooth pairs in contact. The generalized 
mesh force components t and W fc are considered here to be known and fixed. 

Equations (2) and (3) are readily solved for the tangential generalized trans- 
mission error component £ t : 



a r - b r 
lb la 


a l b 2 ' b l a 2 


(5.6) 


Let y g j denote either endpoint position y^j or y B j of a generic tooth pair j in 
contact. In order to show that 5 is stationary with respect to small independent 
variations in the endpoint positions y e j , it is necessary to show that (3? t /3y e j) “ 
0 at the true values of the endpoint positions for all such endpoint positions y e j* 
The generalized mesh force components t and W t are to be held fixed while carrying 
out these partial derivatives. 


We shall let primes denote partial derivatives with respect to a single generic 
endpoint position y e j* Forming the partial derivative of equation ( 6 ) with respect 
to y e j gives 


3 ‘, 

8y 


SYV^VMV b i r i" b i r a ) 
: ( Vb~ b iV <a iY a i b 2 ~ b iY b i a 2 ) . 


ej 


(a iV b l a 2 ) 


(5.7a) 


‘VMWi-NV ' ! t <a i b 2 +a i b 2' b l a 2‘ b ; a 2 > 


a i b 2 ' b i a 2 


(5.7b) 


a i ( Y b 2 t t ) ‘ "I'VVt 1 ' VYYt' + a l (r b‘ b 2 C t ) 


a i b 2 ’ b i a 2 


(5.7c) 
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where equation (6) was used in going from equation (7a) to equation (7b). However, 

from equations (2) and (3) we have r -at = a y and r - b t = b u. 

a 2 t 1 b 2 s, t 1 ’ 

respectively, which when combined with equation (7c) gives 


_ V r ;- b 2V W 

ay 


ej 


a b -b a 
12 12 


(5.8) 


From equations (3.37), (3.38), (3.36a) through (3.36o), and (3.32a) through 
(3.32o), it follows that the coefficients and right-hand sides of equations (2) and 
(3) can be expressed as 


" I 


lj ’ 


l 


V » 

2j 


r = x 
a 


I n 


aj 


(5.9a,b,c) 


and 


b > ‘ \ b lJ - b 2 - J b 2J • r b ' “t + ? "bj > (5.10a, b,c) 

J j 1 

where equations (4) and (5) have been used in equations (9c) and (10c). Each term 

an, a?.,, n ., b tJ , b„., and n. . in equations (9a,b,c) and (10a, b,c) is dependent on 

AJ aj lj 2j bj 

only the endjoints y^. and ygj of the single mating tooth pair j. Moreover, these 
six terms are independent of x and W t « Hence, continuing to denote partial 
derivatives with respect to y g j by primes, we have by differentiating equations 
(9a,b,c) and (10a, b,c) with respect to y e j while holding x and W t fixed. 


a = a , 

1 lj 


2j 


(3r /3y .) = n* . 

a ej x aj 


(5.11a,b,c) 


and 


b' = b' , 

1 lj 


b' = b' , 
2 2j * 


( 3r b /3i, ej )w t ■ "w > 


(5. 12a ,b ,c) 


where the subscripts x and W t denote that these quantities have been held constant 
while evaluating the derivatives. Combining equations (11a, b,c) and (12a, b,c) with 
equation (8) gives 


3C. 


lKj~ b 2j C t" b lj y) " V n lj a 2j C t a Ij y ) 


ay 


ej x ,W. 


a b -b a 
12 12 


(5.13) 
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where the subscript x,W t on the left-hand side denotes that j and W t have been held 
constant while evaluating the partial derivative ^^ t ^y e y 

From equations (3.12) and (3.10), we see that the generalized mesh force 
components x and W t can be expressed as 


l 

j 


j 


w = l 

t L 

j 


w 


tj 


(5. 14a ,b) 


where Xj and W t j are the components of x and W t , respectively, transmitted by tooth 
pair j. It follows from the discussion between equations (3.33) and (3.37) and 
equations ( 2 ), (3), (9a,b,c), (10a, b,c), and (14a, b) above that Xj and 
W t j can be expressed as 


a u 

lj 


+ a r — q 
2j t aj 


(5.15) 


and 


W 

tj 


= b y 

lj 


b ? 
2 j^t 



(5.16) 


which are abbreviated forms of equations (3.33) and (3.34), respectively. Forming 
the partial derivatives of equations (15) and (16) with respect ot y e j while holding 
the generalized transmission error components p and 5 constant gives 


and 



V >5. 


a' u 

lj 


+ a' r 
2 j c 



(5.17) 


3W 


3y 


ej y,C 


- b > + b :.? 




2 j t 


- n ' 


bj 


where primes denote partial derivatives with respect to y e j , 
equations (17) and (18) into equation (13), we have 


as before. 


(5.18) 

Substituting 


9y ./ TI 

ej ' x ,W 


= - (a b -b a ) 1 
12 12 


(aw 


a 


tj 


ax 


i\ay ./ 

X ej' p ,r 


- b 


j 


lV 3 y aj/p, ?t j 


(5.19) 
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From equations (3.13) and (3.16), respectively, it follows that W t j and Tj can 
be expressed as 


and 


y Bj 

Vf “ J g, 4 (x,y)u (x,y)dy 
tj lj nj 

y AJ 


BJ 

r o / [g 2 j (x,y) - g 1 j(x,y)]u^(x,y)dy , 

y Aj 


(5.20) 


where 


( 5 . 21 ) 


g (x,y) = sea|) (y)K ( ^(x,y)cos[i|» (y) - X.(x,y)] 
lj b nj b j 


(5.22) 


and 


g (x,y) ± [ l-(y/r ^K^U.y) 
2j o nj 


(5.23) 


where r and y are defined in Fig. 3, x is defined by equation (2.6) and Fig. 3, 

o (t) (u ) 

angles ij),(y) and X.(x,y) are defined in Fig. 9, K (x,y) and K (x,y) are local 
b j nj nj 

stiffnesses of tooth pair j per unit length of line of contact defined in Appendix B, 

and u ,(x,y) is the local combined elastic deformation of the running surfaces of 

nj 

both teeth of tooth pair j evaluated at location y along the line of contact, as 
described by equation (3.1) and the accompanying text. 


According to equation (19), y and £ must be held constant when evaluating the 
partial derivatives of the right-hand sides of equations (20) and (21). From equa- 
tion (3.28), it follows that the elastic deformation u n j(x,y) must be treated as a 
constant when y and £ are held fixed. Furthermore, in the steepest descent 
procedure used to compute the endpoint positions y^j and ygj and the generalized 
transmission error components y, and £^, the local tooth pair stiffnesses 
K^^(x,y) and K^j^(x,y) will be assumed to be independent of y^j and ygj during any 
descent. Forming the partial derivatives of equations (20) and (21) with respect to 
the general endpoint variable y e j and treating K^\x,y) and K^\x,y) as independent 
of y^j and y B j » we thus find 
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(5.24) 


3W 


,ay 


tA 

ey y,?. 


± g,,( x »y ,)u ,( x »y ,> 
lj ej nj ej 


and 


3t 


3y 


ej u,t 


* t ot g 2 j U - y eJ ) ' g l 3 (x - y ej ) J U „ j U - y eJ ) ’ 


(5.25) 


where the upper (plus) signs apply to the case y e j * y B j and the lower (minus) signs 
apply to the case y e j “ y^j* 

Finally, substituting equations (24) and (25) into equation (19), we obtain 


*1 

3y 


a i g lj" r o b l (g 2 j' g lj ) 
= + 


ej t,VJ 


V2-V2 


y=y 




(5.26) 


ej 


where the upper (minus) sign applies to the case y e j ■ y B j and the lower (plus) sign 
applies to the case y e j « yA j * The quantities a^, a^, b^, and b 2 are the co- 
efficients in equations (2) and (3) which are readily determined by comparing these 
two equations with equations (3.37) and (3.38), respectively. The quantities gjj and 
g 2 j are defined by equations (22) and (23). Each of these six quantities is 
dependent on the endpoint value y g j under consideration and is to be evaluated at 
that value in equation (26) as is indicated. 

For a fixed rotational position of the meshing bevel gear pair designated by the 
variable x, as we proceed along a generic line of tooth contact j from the interior 
of the contact region outward in either possible direction, the composite local 
elastic deformation u n j(x,y) of the tooth pair diminishes from a positive value in 
the interior of the contact region to a value of zero at the boundary y = y g j of that 
region. Thus, when the local tooth pair stiffnesses K^j^(x,y) and K^^(x,y) are 
considered to be independent of the endpoint positions y = y A j and y = y B j of the 
line of contact of tooth pair j, equation (26) shows that the tangential component 5 
of the generalized transmission error is stationary with respect to small variations 
in the endpoint positions of the line of contact of tooth pair j about their true 
values when the components t and Wj. of the generalized force transmitted by the mesh 
are held constant while the endpoint positions are varied. Furthermore, since 
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equation (26) applies to any tooth pair j in contact, we have shown that g is 
stationary (ref. 22) with respect to small independent variations in the line of 
contact endpoint positions of all tooth pairs in contact. 


Stationary Property of y(x). The symmetry in the relationships leading to 

equation (26) enables us to write an expression for (3y/9y ,') directly from the 

e j t ,w 

result (26). To see this, we note first that the derivation of equation (19) is 
based entirely on the relationships described by equations (2) through (5), (9a,b,c), 
(10a, b,c), (14a, b), (15), and (16). If the members in each of the five pairs of 
quantities (y,g^), (a^,a^), (b^,b^), (a * a 2 j ^ * an< * ^lj » ^2j ^ are exc ^ an S et * with one 
another, the above set of relationships remains unchanged except for some of the 
equation numbers. Making these same exchanges in equation (19) yields 


3y 


3y 


ej / t ,W 


- ^VYi^ 1 


2 W ej A t >y 2 \ 3y ej/ c t ,y 


(5.27) 


where the order of the quantities and y in the subscript £ t ,y is of no conse- 
quence. Furthermore, equations (20) through (25) remain unchanged when the members 
in each of the above-mentioned five pairs of quantities are exchanged with one 
another. Hence, we obtain a valid relationship by exchanging the members in each of 
the above five pairs of quantities in equation (26), which yields 




a 2 g l j - r o b 2 (g 2j- g l j ) 


a 2YYl 


y=y 


V x -V • 


ej 


(5.28) 


where the upper (minus) sign applies to the case y e j = y B j and the lower (plus) sign 
applies to the case y e j = y^ j . It follows from equation (28) that the rotational 
component y of the generalized transmission error is stationary with respect to small 
independent variations in the endpoint positions of the lines of contact of all tooth 
pairs in contact about their true values for the same reason given in the discussion 
of equation (26). 


Stationary Property of C r (x). Equation (3.39) is of the form 


cy + c r = W + n , 
1 3 r r c 


(5.29) 


70 


where £ r is the third component of the generalized transmission error and W r is the 

third component of the generalized force transmitted by the mesh. Expressions for 

c , c , and p are obtained by comparing equation (29) with equation (3.39). Each of 
X 3 c 

these three terms is an explicit function of the endpoint positions y A j and y B j of 
all of the tooth pairs j in contact, as one can see from equations (3.39), (3.36a) 
through (3.36o), and (3.32a) through (3.32o). Thus, y and 5 ^ are functions of these 
endpoint positions. Furthermore, from equations (2) through (5) we have seen that y 
must be regarded as an implicit function of the generalized mesh force components 
t and W^. It follows then from equation (29) that must be regarded as a function 
of all three mesh force components t , W , and W^. 


Forming the partial derivative of equation (29) with respect to a single generic 
endpoint position y e j while holding all three mesh force components t, W fc , and W r 
constant, we have 




i lay 


ej/ t ,w 


+ cjy + c 3 ( 


ej ' t ,W ,W 
t r 


+ c's 

3 r 


n . 

c 


where primes denote partial derivatives with respect to y e j , as before, 
the terms in equation (30), we have 


(5.30) 

Rearranging 


(— — \ - ~ c l— — \ “ (c’y + c'c -n') 

\3y / 1\ 3 y / ^ 1 3 r c J 

V ej / t,W ,W \ ej /t,W 
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ej/ t » w 
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(5.31a) 


(5.31b) 


ej' u,C. 


t " 'r 

where equation (29) was used in going to the second line and subscripts following 
parentheses in equations (30) and (31a,b) denote quantities to be treated as 
constants in evaluating the partial derivatives. 


From equation (3.11), we see that the component W r of the generalized mesh force 
can be expressed as 


W 


l W . . 

j rJ 


(5.32) 


Furthermore, from equation (3.35) it follows that W r j can be expressed as 


W 

rj 


c y 

u 


+ c r 
3j r 



(5.33) 
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where by comparing equations (3.35), (3.36a) through (3.36o), and (3.39) with 
equations (33), (32), and (29) above, we see that 


c - Y c , 

1 1 u 


c = y c , 

3 * 3J 


n ■ l n » 

c J CJ 


(5.34a,b,c) 


where each term , cgj , and n c j is dependent on only the endpoint positions y^j and 
y B j of the single tooth pair j. Hence, from equations (29) and (32) through (34c) it 
follows that 


aw 


aw 


rj 


dy . \ay , / 

ej 7 y ej ' v»? 

r r 


(5.35) 


From equation (3.14), it follows that W r j can be expressed as 


'Bj 

W . = / g (x,y)u (x,y)dy , 

rj J 3 nj 

Aj 


where 


,(r) 


g.(x,y) = sectj; (y)K . (x,y)sin[i|> (y) - X (x,y)] 
3 b nj b i 


(5.36) 


(5.37) 


where the local tooth pair stiffness K , \x,y) is defined in Appendix B. 

n J 

Furthermore, it follows from equation (3.27) that u n j(x,y) is independent of y^j and 
ygj when y and c r are held constant. Hence, considering K^^(x,y) to be independent 
of y^j and yg j > 1* follows from equations (36) and (37) that 


'aw 


rj 
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ej ' v 
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3 ej nj ej 


(5.38) 


where the upper (plus) sign applies to the case y g j ■« y B j and the lower (minus) sign 
applies to the case y e j = y^j* 

Finally, combining equations (31b), (28), (35), and (38), we have 


(M .*j_l 

1. 

a 2®lj~ r o b 2^®21 _ ®lj > " 
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where the upper (plus) sign applies to the case y e j = y B j and the lower (minus) sign 
applies to the case y e j = y A j. It follows from equation (39) that the radial 
component ? r of the generalized transmission error is stationary with respect to 
small independent variations in the endpoint positions of the lines of tooth contact 
about their true values for the same reason given in the discussion of equation (26). 


Minimum Property of Generalized Transmission Error Component 

In order to show that the tangential component 5 ^ of the generalized 

transmission error takes on a minimum value when the endpoint positions y A j and y B j 

of the lines of contact of all tooth pairs j in contact take on their true values , 

it is necessary to consider the second variation (ref. 22 ) in r at the true values 

t 

of the endpoint positions. For 5 ^ to take on such a minimum it is necessary that 
the second variation be positive, which requires (ref. 22 ) in our case that 





3 V y ek/ x ,W 


a a 
i k 


> 0 


(5.40) 


for all real numbers and subject to the condition for all j that | ot ^ j < 1 and 
that Z <Xj 2 =1. In interpreting y e j (or y e ^) here, endpoints y A j and y B j on the 
same 4ooth pair are treated as y e j for two different values of j; i.e., j and k 
count endpoints here, two per tooth pair, not tooth pairs. The condition (40) 
clearly is satisfied if 


* 2 *t 


^ej 9 Wx,W t 


k * j 


(5.41) 


and 


3 2 r 




> 0 


(5.42) 


T >W 


for both endpoint values y e j = y A j and y e j = y B j of all tooth pairs j in contact. 
It is shown below that the condition (41) always is satisfied, and the condition 
(42) usually is satisfied and always can be satisfied by a readjustment of the 
parameter r Q . 
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Mixed Partial Derivatives of 5 ^. Equation (13) is of the form 
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= - = D-lN . 
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Hence , 


8 V y e k/,,w 


2D 


ay 


ek 


9y 


ek 


(5.43) 


(5.44a) 
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9y ek 9y ej 9y e k> 


(5.44b) 


where equation (43) was used in proceeding from equation (44a) to equation (44b), 

and t and are to be treated as constants in evaluating the partial derivatives on 

the right-hand sides of these equations. Since from equation (26) it follows that 

/3y ,) = 0 at the true endpoint positions, it is necessary here to show only 

t ej t ,w 

that [9N/3y ek J T ^ = 0, provided that D = a^b^ - b^a^ * 0. 

Continuing to leave off the explicit notation that t and W t are to be held 
constant in forming derivatives, it follows by differentiating the numerator in the 
right-hand side of equation (13) with respect to y ek that 
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(5.45) 
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where we have recognized that the primes in equation (13) denote partial derivatives 
with respect to y ... Since the vanishing of 3r /3y . and 3p/3y . at the true 

e J t 6J SJ 

endpoints implied by equations (26) and (28), respectively, applies to any endpoint, 
we have 



3y _ 


3y 


= 0 


ek 


(5.46a,b) 


at the true endpoints. Furthermore, from the vanishing of u n j(x,y e j) at the true 
endpoints it follows from equations (17) and (25), and (18) and (24) that at the 
true endpoints we have 



-at 
2j t 



v = o. 


bj 
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2j t 
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(5.47a,b) 


Thus, for 3N/3y gk to vanish for different endpoints y e j and y e k> it further is 
required only that (a 2 ^ /3y ej 3y ek ) , ( 3 2 b 2j /3 y g j 3 y ek ) , (S^^/Sy^Sy^) , 

(3 2 n aj / 3 y e j ayek)» ( 92a 2j /9 y e j 3y ek) » and ^ 2a lj /3y ej 3y ek) each vanish. Comparing 
equations (2), (3), (9a,b,c), and (10a, b,c) with equations (3.37), (3.38), (3.36a) 
through (3.36o), and (3.32a) through (3.32o), we see that each termri b j, b^j , bj j , 
r\ aj, ^ 2 j , and ajj is a linear combination of integrals of the general form 

A ^ 

H ^ X ’ y Ai’ y B^ = f h(x,y)dy , (5.48) 

d y 
Aj 

where h(x,y) is independent of all endpoints y^j and y B j when K^j^(x,y) and K^j^(x,y) 
are regarded as independent of y^j and y B j. Hence, each term » ^2 j » bjj , h a j » 
a 2 j * and a^j also is of the form of equation (48). However, differentiating 
equation (48) with respect to either y e j = y^j or y e j * y B j , we obtain 


3H 


3y 


ej 


ih(x,y e j ) ; 


hence, whenever y ek * y g j , it follows from equation (49) that 


(5.49) 


3 2 H 

8y ej 9y ek 


= 0 


y ek * 


'ej 


(5.50) 


Therefore, whenever y ek * y e j * we have shown that (3N/3y ek ) = 0 at the true 
endpoints, and from this result and the vanishing of 9 C t / 9 y e j at tlie true endpoints 
it follows from equation (44b) that 
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I = 0 , y * y 

5 VWt,W 6lC 6J 
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Equation (51) is valid at the true endpoints of the lines of contact, provided that 
D = a^b 2 “ bj^ * 


at the true 


Unmixed Partial Derivative* Of e t - If we set y ek = y e j in equation (44), we 

can use that result to evaluate (a 2 £ /3y . 2 ) n . Since ( 3 ^ / 3 y ) = 0 at the trt 

t ej t ,w t ^ t ey 

endpoints, we proceed directly to the expression (45) with y ek = y e j • Again 
utilizing equations (46a, b) with y ek « y e j and equations (47a, b), we find at the 
true endpoints that 


“iKj - b 2j C c ' - b lKj ‘ a 2J C t ' a i' j u) ■ 


-{Pi) 

\ e J / \ ej / u ,5 


where equations (18) and (17) were used in proceeding to the second line. The 
right-hand side of equation (52a) is an explicit function of neither T nor 
Differentiating equations (24) and (25) with respect to y^ and evaluating the 
resulting expressions at the true endpoints where u * = 0 , we obtain 


W7) ' “ll"* ’ 

ej / 


y 1 = (g -g )u' . , 

3y / / o 2 j lj nj 

ej /P,? 


where the primes again denote differentiation with respect to y e j* ^et us define 


A.(x,y .) 
2 ej 


.(x,y .) - r b [g (x,y .)-g,.(x,y .)] 
1 ej o 1 L 2 j ■'ej °lj e.i J 


A 1 lj ej 
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Then, substituting equations (53a, b) into equation (52b) and the resulting expres- 
sion into equation (44b), and recognizing that ( 3 £ / 3 y .) = 0 at the true endpoints 
and that D = a^b^-b^a^, we us ^ n S the definition (54) that 
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(5.55) 
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where the upper (minus) sign applies to and the lower (plus) sign applies 

to y e j = y^j . Equation (55) is applicable only at the true endpoints y^j and y B j 
where the elastic deformations u n j(x,y e j) vanish . The result given by equations 
(54) and (55) also can be obtained formally by differentiating equation (26) with 
respect to y g j and setting u n j ■ 0 in the resulting expression. 


For to take on a minimum value at the true endpoints y^j and y B j> equation 
(42) must be satisfied at each true endpoint position. This criterion is readily 
understood if we consider the Taylor's series representation of as a function of 
endpoint position y e j about the true endpoint position y°j» while holding all other 


endpoint positions constant, 
this representation is 


Retaining terms through the second degree in y e j “ y^ , 


- r 

e ■ 5. n + c'(y-.-y°) + 7 c" (y -y ° .) 

t tO tO ej ej 2 tO ej ej 


(5.56a) 


1 2 
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So 2 So ej ej 


(5.56b) 


where primes denote derivatives with respect to y g j while holding x and W t fixed, 

and subscripts 0 denote that 5 and its derivatives are evaluated at the true 

endpoint position y° . Since from equation (26) we have r’ = fgr /ay ) =0 

ej t ^ t ej' , x,W t 

at the true endpoints, the term proportional to 5 ^ is missing from equation (56b). 

For to increase as y e j is varied in either direction about its true value 

y 0 ^, it is clear from equation (56b) that 5 ^ must be positive. Hence, the 

condition (42) requires that the right-hand side of equation (55) be positive for 

every endpoint. 


This requirement is readily understood with the aid of Fig. 16, which shows the 
behavior of the elastic deformation u^tx.y) of tooth pair j near the true endpoints 
y^j and y B j of the line of contact. This elastic deformation is denoted by the 
solid lines. Let us assume, now, that the local stiffnesses K^5^(x,y), K^^(x,y), 


,<r). 


nj 


nj 


and K" (x,y) in equations (37), (23), and (37) are extrapolated beyond the true 


contact region y^j < y < y B j using the stiffness values K^ 3 (x,y a .. ) , 
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( r ) 

and K nj ( x »y ej ) at the endpoints, y g j = y A j or y B j , as appropriate. Then, if the 
endpoints y A j or y B j used in equations (3.37) through (3.39) are chosen outside the 
true contact region, the equations will predict the dashed behavior of u n j(x,y) 
illustrated in Fig. 16, resulting in a prediction of tensile forces between the 
teeth outside the true contact region. This predicted fictitious behavior can be 
put to good use in determining the true endpoint positions. 

Let us now use this fictitious extrapolated behavior of u n j(x,y) in the 

neighborhoods of y = y A j and y = y B j in our interpretation of equation (55). From 

Fig. 16 at y ■ y B j we have u^ < 0, whereas at y = y A j we have u^j > 0. Utilizing 

these signs of u ^j > it follows from the sign convention of equation (55) that 

( 9 2 £,./ 3 y . 2 ) „ will be positive if 

l ej x 

Aj(x,y g j) > 0 (5.57) 

at all true endpoint positions y e j = y A j and y e j = y B j , where we recall that 
Aj(x,y e j) is defined by equation (54). 

The requirement (57) therefore replaces the requirement (42) for 5 to take on 
a minimum. It is shown in Appendix F that the requirement (57) usually will be 
satisfied, and that by appropriately adjusting the value of r Q the requirement ( 57 ) 
always can be satisfied. In the computational procedure described next, it is 
assumed that Aj(x,y e j) > 0 since this requirement always can be satisfied by an 
adjustment in the value of r 0 « 
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FIG. 16. PREDICTED ACTUAL (SOLID LINES) AND FICTITIOUS (DASHED LINES) LOCAL ELASTIC 
DEFORMATION OF TOOTH PAIR j NEAR ENDPOINTS y Aj AND y Bj OF LINE CONTACT OF 
TOOTH PAIR j. 
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COMPUTATIONAL PROCEDURES 


A method is outlined below for numerical computation of the endpoint locations 
y^j and y B j of the lines of tooth contact and the generalized static transmission 
error components ? t (x) and y(x) at a fixed rotational position of the gear pair 
described by x = x k . The method is based on the intuitively appealing method of 
steepest descent (ref. 25) sometimes called the gradient method. 


Line of Contact Endpoints and Generalized Transmission Error 
Components for a Fixed Rotational Position of the Gear Pair 


The method is based on a very simple idea. For a given rotational position of 
the gear pair described by the value x * x k , the generalized transmission error 
component ^(x^) considered to be a function of the two endpoint vectors 
{y Aj Cx^)} and {y B j( x k )}> where the number of elements j in each vector is the number 
m of tooth pairs j in contact. Thus, we can envision ^(x^) to he displayed in a 
space of 2m+l dimensions as a function of the 2m variables y A j( x k)» yBj^ x k^» j " jp 
^2 » • • • »Jm* 


Trial vectors {y Aj ( x k >} and {y Bj ( x k )} are chosen. A suitable choice for these 

trial vectors is the endpoint values y A j( x k) and yBj^ x k^» J * ji»j2»***»3m 
determined by the geometric zone of tooth contact illustrated in Figs. 3 or 7, as 
appropriate. Since it was shown above that ^(x^) takes on a minimum value when the 
components of fy^^ 3 ^)} ant * ^Bj^ X k^ are 8* ven their true values, a new trial set 
of endpoint values is computed next , where the new set is obtained by proceeding in 
the direction in the 2m+l dimensional space that maximizes the rate of descent 
of ? t ( x k ) from its initial trial position. This direction is the negative of the 
gradient of ^(x^). 


-VC t - (-V 3y A J - ' - V 3 Vi ' 55 t /3y Bj 

1 m 1 m 


( 6 . 1 ) 


The components 3s t /3y A j and 3? t /3y B j, J “ are computed by equation (5.26), 

where u n j(x,y) is determined by equation (3.28) after employing equations (3.37) and 
(3.38) to compute ^(x^) and l^ 3 ^) use * n equation (3.28). A new pair of 
generalized mesh force components W t (x k ) and x(x k ) is to be used in the determina- 
tion of ^(x^) an d u(x^) by equations (3.37) and (3.38) for each new set of 
components an< * 3? t /3y B j , j = as discussed in the next subsection. 
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A set of 2m increments for the m values of j is computed next. 



-h9r /9y , 

t Aj 
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where for all values of j 

increments 6 y.. and 6 y . 
Aj 

( yBj-yAj>* J ■ Jl Jm* 


“ h is a constant chosen so as to keep 

small in comparison with the largest of the m values of 
New trial endpoint vectors are then determined by 
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(6.3b) 


where i counts iterations. These values provide a new set of trial vectors 
{y^ + ^(x^)} and ^ rom which new values and y^ + ^(x^) 

are computed by equations (3.37) and (3.38) using the same generalized mesh force 
components that were used in computing the derivatives 9? t /9y A j a °d 3 C t /dy]$j in 
equations ( 2 a) and ( 2 b), respectively. 


The new value ^(x, ) should be smaller than its predecessor 5 ^ ^(x, ), since 

t x / ... . / J \ c K 

the negative of the gradient points downhill. If 5 ^ '(x^) < '(x^), the above 

procedure is repeated again with the same value of h, starting with recomputation of 
the gradient components 9 £ t / 9 y A j , 9 c t / 9 y B j, J “ ji,...,j m us * n §> new values of the 
generalized mesh force components, as described above. This computation will result 
in a new value 5 ^ + ^)(x^), and so on. As long as c^^^x, ) < ^(i+^~D ( x ^) t the 
same value of h is used to produce each new value of However, at some 

iteration the newest value of ^(i + ^)( x ) w m be larger than the prior value 

/ J I Q _ 1 \ t K 

5 ^ ” '(x^). This will happen when the steps 6 y A j and 6 y R j» j “ jj,...,j m computed 

by equations ( 2 a) and ( 2 b) overshoot the minimum value of ^(x^) i- n our 2 nH-l 
dimensional space. 


(i>< 


When the new value of C^ + ^^(x^) exceeds 


substantially reduced, and new values of Sy 


Aj 


^ (i+ 8,-1 )( x \ ^ t he va i ue Q f b is 
t K 

and 6 y Bj , j = j x , . . . , j m are recomputed 
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by equations (2a) and (2b) using the last set of already existing derivatives 3? t /3y^j 

and 3 £ /Sy,,.* j = . The iteration procedure then is continued as described 

t Bj 1 m . . 

above until an increase in the value of j" ^'(x^) again is obtained. At this 
juncture, the value of h again is substantially reduced. 

This procedure is continued until the absolute value of the largest of the set 
of 2m increments 6y A j and 6y B j , j = that produces an increase in 

^^(x^) is negligible (by several orders of magnitude) in comparison with a typical 
value of (y B j - y A j). When this event occurs, the minimum value of ^(x^) in our 
2m+l dimensional space effectively has been reached. 


Since both $ t (x) and p(x) are required to compute u n j(x,y) using equation 

(3.28), both £ t (x^) an d iKx^) are available at the termination of the above 

procedure. The values hCx^)* an d the final endpoint vectors fy^j(x^)} an< * 

{ysj( x k)}> i = are t0 stored. The third generalized transmission error 

component £ (x,) can be computed from the final endpoint values and p(x. ) using 

equation (3.39). The resulting triad of transmission error components p(x^), 

£ (x ), and r (x. ) is the generalized static transmission error evaluated at the 

t k r k 

rotational position of the gear pair described by x = x^. 


When the above procedure is started, care must be taken to insure that the set 

j = jj j m includes all tooth pairs in potential contact. At each iteration, 

however, there is no guarantee that the new set of endpoints computed by equations 
(3a) and (3b) will satisfy 


(i+1) > (i+1) 

y Bj ' y Aj 


j " J i»* ,, » J m 


(6.4) 


For each value of j not satisfying equation (4), one must set = yfV"^ in the 
remainder of that iteration, thereby eliminating those lines of tooth contact in the 
computation of 5^ + ^(x^) and y^ + ^(x^). 


In related problems, we have found the above-described gradient descent to work 
very well. It requires a minimum computer programming effort. 


Generalized Mesh Force Components 


the 


For the same rotational position of the gear pair described by x = x^, 
generalized mesh force components W^^x^) and are rec l u: *- re< l f° r computation 
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of and y^(x^) by equations (3.37) and (3.38) using the endpoint vectors 

{ y Aj (x k } J and ^ y Bj )(x k^’ j = as described above. These generalized 

force components are to be computed using equations (4.24) and (4.19). However, the 
coefficient matrices K(x) and n R (x) in equation (4.24), which also appear in the 
matrix form (3.40) of equations (3.37) through (3.39), depend on the endpoint 
vectors {yj^(x)}, {y^C*)}* j = j i* - • - . J m - Hence, generalized mesh force 

components w£^(x^) and T^ d ^(x^) computed using the endpoint vectors {y^^x^)} anc * 

{yg^(x^)} should be employed in the above-described gradient descent procedure. 


The sequence of steps involved in the iterative portion of the procedure for 
computing the line of contact endpoints and generalized transmission error 
components for a given rotational position x = x^ of the gear pair is listed below: 


A. Initial trial endpoint set {y A j^}> {y B j^}» ^ = ^l , ***’^m * 

B. Compute and x^ with equations (4.24) and (4.19) using line of 

contact endpoints {yj^} , {y^^}. 

C. Compute , y^, and u^^ at endpoints with equations (3.37), (3.38), 
and (3.28) using endpoints {y^} , { yjg j ^ • 


D. 


E. 


Compute 3C t /9y A j and 3S t /3y B j> j - J i» - - - . J m with equation (5.26) using 
endpoints {yj^}, {yj^}* 

Compute new endpoint values {y^ d+ ^}, {y B j + ^}> j " jj j m with 

equations (2a,b) and (3a,b) above. 


F. 


G. 


Test: Is y£j +1) > y A j +1) . J = 5 v 




If not, set y 


(i+1) - „(i+D 


Bj 7 Aj 

Compute and y^ i+1 ^ with equations (3.37) and (3.38) using new 

endpoint values {y A j + ^} , {y B j + ^} » j = and values of 


W^ i) and x^ d ^ computed in step B. 
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I 


H. Test: Is • If yes, delete increase index i 

by one and proceed to step B. If no, delete c£ i+ ^» and test 

maximum increment magnitude 1 6 j ^ | or scanne< * over all j = 

j 1» • •• , j m » where 6yj^ and <$Ygj\ 3 = J i » • • • » J m are those values used in 

step E. If maximum increment magnitude is negligible (several orders of 

magnitude) in comparison with typical interval (yif**^ - y$* + ^), effective 

A3 


minimum value of C t “ has been achieved. If maximum increment 


magnitude is not negligible, substantially reduce h and proceed to step E« 

Comment: Notice that the comparison test "Is always is made with 

(i+1) (i) C t 

values of and 5^ having been computed using the same pair of generalized 

C (i) (i) 

mesh force components y and x since minimum property of true value of 5 was 
proved under assumption of constant mesh force components W t and x. 


Generalized Transmission Error Components and their Spectra 


The above procedure generates the generalized static transmission error 
components y(x), C t (x), and ? r (x) at a single rotational position of the gear pair 
described by x = x^. By repeating this procedure at an equispaced sequence of 
values of x^> the generalized transmission error components (and lines of tooth 
contact endpoints) are determined as functions of the rotational position of the 
gears. 

If the elastic properties and running surface modifications from perfect 
equispaced spherical involute surfaces are the same for every tooth on each gear of 
the gear pair, then the generalized static transmission error components y(x), 

£ t (x), and C r ( x ) each is a periodic function of x with period A equal to the tooth 
spacing interval in x (Fig. 3 ). Thus, minimum definitions of y(x), ? t (x), and ^(x) 
as functions of x would require the above described procedure to be caried out at 
about 8 points with spacing in x of A/8. An optimum method of interpolating between 
such equispaced transmission error component computations is described below. 

A representation of the transmission error in the frequency domain (refs. 7 - 9 ) 
often is of primary interest. Such a representation can be obtained from the trans- 
mission error computed at equispaced intervals using the method of trigonometric 
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interpolation (ref. 31). Let us denote a generic component y(x), ? t (x), or 5 r (x) 
by e(x). Assume that ?(x) is periodic with period A, as noted above, and 
that s(x) has been computed at 2n equispaced points x^, where the spacing between 
adjacent points is 



Let the 2n points x^ at which ^Cx^) has been computed be numbered k = -n, 

-(n-1) , . . . ,0, . . . ,n-l. Since ?(x) is periodic with period A, the value of a 2n+l st 
point k = n is provided by the point at k * -n; i.e.. 


x = x . (6.6) 

n -n 

Let x* denote the independent variable in Lanczos' (ref. 31) treatment of 
trigonometric interpolation. The period of the periodic function there has been 
normalized to the interval [ — ir ,ir ] . To change the interval [— ir,ir] to the period A 
in x of our independent variable, we require the transformation 


x' = 2irx/A . 


(6.7) 


Incorporating equation (7) into Lanczos' (ref. 31) trigonometric interpolation 
formula (4-11.24) gives 


n 


C(x) = £ ' a exp(12irpx/A) , 


( 6 . 8 ) 


P=~n 

where the expansion coefficients ctp are provided by Lanczos' formula (4-11.25), 
n , 


1 ’ 

<* = ^ 4 ?(x k )exp(-i2Trpx k /A) , 

v k=-n 


(6.9) 


where obvious changes in notation have been used. The prime in equation (8) 
designates (ref. 31) that the two limiting terms a_ n and a n are to be weighted by 
1/2 in the summation; similarly, the prime in equation (9) designates that 
?(x_ n ) and ?(x n ) are to be weighted by 1/2 in that summation. 

The prime in equation (8) is absent from the corresponding equation (4-11.24) 
in reference 31. However, the validity of the interpolation formula (8) readily can 


85 


be established by substituting equation (9) into equation (8), interchanging the 
order of summation, and then utilizing equation (4-11.17) of reference 31. 

Particular attention must be paid to the case where equation (8) is used at x = x_ n 
or x = x n » It follows that the interpolation formula is valid at these endpoints 
only if £(x_ n ) = C( x n )> which is satisfied in our application since ^(x) is periodic 
with period A. By substituting Lanczos' equation (4-11.20) into his equation (4- 
11.22), one can see that Lanczos* equation (4-11.22) is not valid for | j — k. | = 2n. 
This same conclusion can be reached by using L'Hopital*s rule on the right-hand side 
of Lanczos* equation (4-11.17). Nevertheless, with the prime included, we have 
shown that equation (8) is an exact interpolation formula when the a are computed 

p 

by equation (9), provided that ?(x ) = r(x ), which is satisfied in our 

— n n 

application. 


Equation (8) is identical to the complex Fourier series representation of the 

transmission error component s(x) (ref. 7, equation (64)) except that (tooth 

meshing) harmonics p only through order n are included in the representation (8). 

Thus, from equation (8), we see that the coefficients a , I p I « 0,l,...,n-l 

P 

generated from computations of c( x ^) at the 2n equispaced points x^, k = -n,-(n- 

1 ) , . . . ,0, . . . ,(n-l) can be directly interpreted as (approximations to) the complex 

Fourier series coefficients of the transmission error component £(x). Because of 

the prime on equation (8), the coefficients a and a generated by equation (9) 

— n n 

(p = ±n) must be divided by 2 to be interpreted as Fourier series coefficients. 


The coefficients a generated by equation (9) can contain contributions from 

high-order harmonics [pj > n "impersonating” the lower-order harmonics |p| < n (ref. 

32, Fig. 8). Such aliasing contamination of the coefficients a can be minimized by 

P 

choosing the number 2n of equispaced points at which 5(x^) is computed to be 
substantially larger than twice the number of accurately desired harmonics |p|. 


86 



Conclusion 


It has been shown that the traditional definition of the one-component trans- 
mission error is unsuitable for accurately characterizing the vibratory excitation 
arising from meshing pairs of spiral bevel gears. A three-component generalized 
transmission error definition was introduced that is made up of two orthogonal 
translational components and one rotational component of relative displacement of 
the gears with respect to their perfect spherical involute counterparts. All 
displacements that contribute to these three transmission error components take 
place in the plane of tooth contact. A transformation of coordinates, equations 
(2.10) through (2.13), was derived that relates the three generalized transmission 
error components to the six degree-of-f reedom deviations in the positions of the 
gear shaft centerlines from the positions occupied by their perfect spherical 
involute counterparts. The methodology is applicable to bevel gear pairs whose 
actual surface of tooth contact is reasonably well approximated by a plane surface 
passing through the point occupied by the apexes of the two pitch cones of the 
mating gears. 

In the absence of friction, the forces transmitted by the meshing teeth can be 
characterized by two orthogonal force components in the plane of tooth contact and 
one torque component, where the vector describing this torque component is normal to 
the plane of contact. Equations were derived that relate the three generalized 
transmission error components to the above-mentioned three generalized mesh force 
components [equations (3.37) through (3.39)]. These equations include contributions 
to the transmission error components arising from elastic deformations of the teeth 
and gear bodies, and contributions from deviations of the running surfaces of the 
teeth from equispaced perfect spherical involute surfaces. A method was described 
for incorporating these equations into the equations of motion of a gear system. 

When the inertial forces of a meshing gear pair and bearing supports are 
ignored, the three generalized mesh force components are determined by the tooth/ 
gearbody stiffnesses, deviations of the tooth running surfaces from equispaced 
perfect spherical involute surfaces, gearshaft bearing support flexibilities, devia- 
tions of the bearing centerlines from the positions occupied by the shaft center- 
lines of rigid perfect spherical involute gears, and the shaft input torque. 
Equations were derived providing the three generalized mesh force components in 
terms of the above-described quantities [equations (4.24) and (4.19)]. The three 
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components of the generalized static transmission error are predicted by equations 
(3.37) through (3.39) or equation (3.44) when the generalized force components 
provided by equations (4.24) and (4.19) are used in their evaluation. 

Finally, a method was developed for computing the positions of the endpoints of 
the lines of tooth contact in the plane of contact. These endpoint positions are 
needed for evaluation of integrals that are required to evaluate coefficients in the 
above-mentioned equations. The computational procedure developed in the last sec- 
tion provides the three components of the generalized static transmission error at a 
set of equispaced points taken over one tooth spacing period of the rotational 
positions of the gears. Equations (6.8) and (6.9) may be employed to provide a 
smooth interpolation between such points for each of the three generalized trans- 
mission error components. When the running surfaces and stiffnesses of all teeth on 
each gear of a meshing pair are the same, equation (6.9) may be used to provide the 
Fourier series coefficients of the tooth meshing harmonics of each of the three 
generalized transmission error components. 
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APPENDIX A 


DERIVATION OF EXPRESSION FOR GENERALIZED TRANSMISSION ERROR IN TERMS OF 
GEAR DISPLACEMENTS AT SHAFT CENTERLINES 

(.) (•) (•) 

We first express the gear body shaft centerline translations Xj » x 2 » x 3 *- n 

the y , y , y coordinates shown in Fig. 6. The y coordinates are obtained by a 
1 2 3 / \ ✓ \ 

rotation about the x, axis through the angle yi* • The common plane of the 
x 2 » x 3 , y 2 » and y 3 axes is shown in Fig. 17(a), from which it follows directly 
that 

*1 " X 1 


y 2 - x 2 cosY b ~ x 3 siny b 
y 3 - x 2 8in Yb + x 3 cosY b , 
or in matrix notion, 



(Al) 


(A2) 


where the first line in equation (Al) is a consequence of the fact that the rotation 
was taken about the x[* ^ axis as indicated in Fig. 6, and where the superscripts 
(• ) have been left off. Next, we express the y 1 , y 2 , y 3 coordinates in terms of the 
z , z , z, coordinates shown in Fig. 6. The z coordinates are obtained by a rotation 
about the y 2 axis through the angle v v * . The common plane of the yj^ , y 3 , , and z 3 

axes is shown in Fig. 17(b), from which it follows that 


z 1 = yjCosv + y 3 sinv 

z 2 = y 2 

z 3 = -y 1 sinv + y 3 cosv , 


(A3) 
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FIG. 17. (a) RESOLVING x 2 AND INTO y 2 AND y 3 COMPONENTS, 

(b) RESOLVING y x AND y 3 INTO Zj AND z 3 COMPONENTS. 
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or 


cosv 


0 


Z 1 

z 2 

Z 3 


! 


0 

-sinv 


1 

0 


sinv 

(M 

0 

r 

cosv 

(y 3 ) 


(A4) 


Substituting equation (A2) into the right-hand side of equation (A4) and multiplying 
the two coefficient matrices yields 



cv 


-sv 


sy sv 
b 

cy sv 
b 

C Y, 

-sy 

b 

b 

sy cv 
b 

cy cv 
b 



(A5) 


where the cosine and sine operators have been abbreviated by c ans s respectively 
(ref. 18, p. 106). 

We are concerned with the four possible gear pair configurations listed in Table 
1. The configuration represented by the first full line in Table 1 is illustrated in 
Fig. 3. When gear (2) in Fig. 3 becomes the driver and its direction of rotation is 
changed from counterclockwise to clockwise, the plane of contact and lines of contact 
shown in Fig. 3 also apply to this configuration, which is represented by the second 
full line in Table 1. Thus, in the configurations described by the first two full 
lines in Table 1, gear (1) in Fig. 3 is represented by the first column in Table 1 
entitled "upper sign" and gear (2) in Fig. 3 is represented by the second column 
entitled "lower sign." Furthermore, the base cone and plane of contact configura- 
tions of both gears shown in Fig. 3 are represented by Fig. 6 which, when upright, 
represents gear (1) and, when upside down, represents gear (2). 

From Fig. 6, we can see that gear body translations in the positive z direction 
are "equivalent" to the addition of material to the running surfaces of perfect 
involute teeth. Hence, for gears (1) and (2), positive contributions to z ^ are 
equivalent to negative contributions to the £ component of the generalized 
transmission error - compare Figs. 5 and 6. It follows from equations (2.10) through 
(2.13) that the last three elements in the middle row of the right-hand side of 
equation (2.13) are the negative values of the three elements in the top row of the 
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coefficient matrix in the right-hand side of equation (A5). Since gear body 
translations in the z 2 direction are normal to the plane of contact, the middle row 
of the coefficient matrix in equation (A5) contributes nothing to the generalized 
transmission error. 

From Figs. 5 and 6, we can see that gear body translations in the z g direction 
provide the contributions to A careful examination of Figs. 3 and 6 shows that 

in the case of gear (1) gear body translations in the positive z direction are 
equivalent to the addition of material to the running surfaces of the teeth on gear 
(1); hence, in the case of gear (1) positive contributions to z g are equivalent to 
negative contributions to the component of the generatized transmission error. It 
follows from this fact and equations (2.10) through (2.13) that in the case of gear 

(1) the last three elements in the third row of the right-hand side of equation 
(2.13) are the negative values of the three elements in the bottom row of the co- 
efficient matrix in equation (A5). Further careful examination of Figs. 3 and 6 
shows that translations of the body of gear (2) in the positive z 3 direction are 
equivalent to the removal of material from the running surfaces of the teeth on gear 

(2) ; hence, in the case of gear (2) positive contributions to z are equivalent to 
positive contributions to 5 ^. Thus, in the case of gear (2) the last three elements 
in the third row of the right-hand side of equation (2.13) are the same as the values 
of the three elements in the bottom row of the coefficient matrix in equation' (A5). 
These facts establish the sign conventions given by the first two full lines in 
Table 1. 

The sign conventions given by the last two full lines in Table 1 are established 
in a similar manner with the aid of Figs. 7 and 8. The counterpart to Fig. 6 in 
these considerations is shown in Fig. 8 which is the mirror image of Fig. 6 - i.e., 
the figure seen by looking through Fig. 6 ftfom its backside. From this fact, it 
follows that equation (A5) also is valid for the gearing configuration shown in Fig. 
7. Consequently, the sign conventions given by the last two full lines in Table 1 
can be derived from equation (A5) and Figs. 7 and 8. Furthermore, from the 
applicability of Fig. 8 to the configurations described by the last two full lines in 
Table 1 it follows from equation (A5) that the last three elements in the second row 
of the right-hand side of equation (2.13) also are applicable to the configurations 
described by the last two full lines in Table 1. 
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We turn now to the four nonzero elements in the upper left-hand corner of the 

matrix in equation (2.13). From equations (2.10) through (2.13), we see that the two 

(.) (.) 

nonzero elements in the first row represent the coefficients of p 2 and p 3 , 
respectively, in their contributions to y(x). Taking into account the sign 


convention for y(x) illustrated in Fig. 5, we can readily see from Figs. 6 and 18 

that the contribution to y(x) from gear (• ) is p^*^ cosy^* ^ - p* ^ siny^ \ from 

1 b 3 b 

which the first row of the matrix in equation (2.13) follows directly. The first two 


nonzero elements in the second row of the matrix in equation (2.13) represent the 


contributions of these rotation components to 5 that arise from the shift of the 
origin through the distance r Q from cone apex to the point 0 ^ on the zone of contact 
centerline shown in Fig. 6 . The signs of these two elements are a consequence of the 


sign convection illustrated in Fig. 5. 


A similar analysis of Fig. 8 , which applies to the gearing configuration shown 
in Fig. 7 and to the last two full rows in Table 1, shows that the four nonzero 
elements in the upper left-hand corner of the matrix in equation (2.13) also apply to 
the configurations described by the last two full rows in Table 1. 
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FIG. 18. RESOLVING (SMALL) ROTATION VECTORS p 2 *^ AND p^ INTO COMPONENTS 

PERPENDICULAR AND PARALLEL TO THE PLANE OF CONTACT. THE PLANE OF THE PAPER 
CONTAINS THE y 2 AND y 3 AXES SHOWN IN FIG. 6. 
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APPENDIX B 


DEFINITIONS AND INTEGRAL EQUATIONS FOR LOCAL TOOTH-PAIR STIFFNESSES K*‘*(x,y) 

n J 

Let kj(v,v';x) denote the composite deformation from both teeth at a point v on 
the line of contact of tooth pair j due to a unit force applied at point v ' on the 
line of contact, where force and deformation both are measured in a direction normal 
to the tooth surfaces and, therefore, in the plane of contact, and x denotes the 
position of the zone of contact as described by equation (2.8). The influence func- 
tion kj(v,v';x) is assumed to include bending, shear, and Hertzian deformation 
components of both meshing teeth, as well as the components arising from deformations 
of the bodies of both meshing gears. The unit force used to define kj(v,v';x) is 
assumed to be applied with uniform intensity to a circular area of diameter equal to 
the width of the contact line between the two gear teeth. Since the Hertzian 
deformation component of real teeth is a (weakly) nonlinear function of loading 
because of dependence of the width of the contact line on loading, specification of 
the area of applied force implies that our influence function is defined for a given 
nominal loading condition. The deformation u^(v;x) at point v caused by a lineal 
force density pj(v';x), where v and v* are on the line of contact, can be expressed 
in terms of kj(v,v';x) by the line integral 


uj(v;x) = / kj(v,v' ;x)pj(v' ;x)dv' , 


(Bl) 


where dv' denotes differential distance measured along the line of contact and where 
the integral is taken over the portion of the line of contact contained in the zone 
of contact. 

Equation (Bl) may be regarded as a Fredholm integral equation of the first kind 
for the lineal force density p^(v';x). We may formally invert (ref. 23, pp. 4-9, 
351-352) equation (Bl) to yield 


Pj(v;x) */ kj 1 (v,v' ;x)u^(v' ;x)dv' , 

where kjMv.v'jx) is the inverse of kj(v,v';x). 

We shall consider K^^(x,y) first. The differential force increments dW n j 
normal to the line of contact are related to the force density Pj(v;x) by 


(B2) 
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(B3) 


dW nj = Pj^ v 5 x ^ dv * 

From equations (3.5), (3.7), and (B3) we can express the contribution from tooth pair 
j to the total tangential force W fc illustrated in Fig. 4 by 

W.(x) = / cos[i|> b (y v )-X.(x,y v )]p (v;x)dv , (B4) 

J v J 

where y^ is an abbreviation for the coordinate y illustrated in Figs. 3 to 5 
expressed as a function of distance v measured along the line of contact of tooth 
pair j. From equation (3.3), we have 

dv = sea|) b (y)dy ; (B5) 

hence, v can be expressed as a function of y by the indefinite integral 

Vy * / seap b (y)dy + C , (B6) 

where C might be chosen by requiring to be zero at y = 0. For a given value of 
zone of contact position x, it is obvious from Figs. 3 to 5 that the function 
Vy * v(y) also determines the inverse function y^ = y(v) which is the function used 
in equation (B4). 

Inserting equation (B2) into equation (B4) and using dv' = sea|> b (y)dy from 
equation (B5), we have 

W t j(x) =/ cos[^ b (y v )-Xj(x,y y )] / kj 1 (v ,v ' ;x)uj (v ' ;x)dv 'dv 


*/ { / cos[i|) b (y v )-X (x.y^jkT 1 (v ,v' ;x)dv} u'(v';x)dv' 

v ' v J J 

y B 

“ / seaj» b (y) { / cos[i|» b (y )-X.(x,y )]kr 1 (v,v';x)dv} 
y A 


x uj (v^;x)dy , 


(B7) 
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where in the last line we have considered v' to be the function of y determined by 
equation (B6). Since UjCv^Jx) and u n j(x,y) both are measured normal to the line of 
contact, we have 


u j (v^;x) = u nj (x,y); 


(B8) 


hence, requiring that equations (3.13) and (B7) represent identical expressions for 
arbitrary deformations u^Cx.y), we must have 

K^(x,y) = {cosU (y)-X (x,y)]} —1 / cos[i|> (y )-\.(x,y )]k“l (v ,v ' ;x)dv , (B9) 

nj b j b v j v j y 

which is the definition of the local tooth-pair stiffness K^^(x,y) per unit length 
of line of contact. The integral with respect to v in equation (B9) is taken over 
the full length of line of contact of tooth pair j at zone of contact position x. 


It is not necessary to determine the inverse influence function k. 1 (v,v ;x) in 

(t) ^ 

order to evaluate the local tooth-pair stiffness Kyy (x,y). The symmetry of the 
influence function k^Cv.v'jx) (Maxwell's reciprocal theorem) implies that kj 1 (v,v';x) 
possesses the same symmetry [23, p. 9] - i.e., 


ky 1 (v ,v' ;x) = ky^v'.vjx) . (BIO) 

Utilizing the property (BIO) in equation (B9) yields 

cos[i|» b (y)-Xj(x,y)]K^^(x,y) = / kyl (Vy,v;x)cos[i|) b (y v )-Xj(x,y v )]dv , (Bll) 

which is of the same general form as equation (B2). Let us define 

u£j(v;x) A cos[ip b (y v )-Xj(x,y v )] (B12) 

and 

P tj (v^;x) A cos[4. b (y)-X j (x,y)]K^ ) (x,y) , (B13) 

which, according to equations (B2) and (Bll), are a deformation/force density pair 
that satisfy equations (Bl) and (B2). It follows, then, from equation (Bl) that the 
solution Pj (v';x) = p t j(Vy5x) to the integral equation (Bl) with left-hand side 
uj(v;x) = u b j(v;x) given by equation (B12) yields the local tooth-pair stiffness 
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(B14) 


K^(x,y) = sec [’*'b (y)_X j (x,y ^ P tj (v y ;x) 


9 


where the function v’ = v'(y) is determined by equation (B6) . 


The situation with regard to the local tooth-pair stiffness K v . (x,y) associated 

nj 

with the radial force component W r j(x) is completely analogous to that above for 
K^j^(x,y). Using equations (3.6), (3.8), (3.14), and (B3), we obtain in this case, 
instead of equation (B9), the definition 


K^(x,y) A | sin[i|> b (y)— Aj(x,y)] 1 sin[> b (y v )-Xj(x,y v )]kjl (v ,v^;x)dv . (B15) 

The local tooth-pair stiffness K^5^(x,y) can be obtained from the solution 
— nj 

Pj(VyJx) = p^(y^;x) of the integral equation (Bl) using for its left-hand side 

u!(v;x) = u' .(v;x) where 
3 rj 


u^.j (v ;x) A sin[i|) b (y v )-X J (x,y v )] 


and 

P r j(Vy;x) - sin[i|» b (y)-X^(x,y)]K^ ) (x,y) 

(r) 

from which one obtains (x,y) _b£ 

K nj )(x,y) = c8c [ , l»b (y) " X j (x * y) ] p rj (v y ;x) * 


(B16) 

(B17) 


(B18) 


where the function v' 

y 


v'(y) is determined by equation (B6) . 


The situation with regard to the local tooth-pair stiffness K^jj^(x,y) appearing 
in equation (3,16) is slightly more complicated. Inserting equation (B3) into 
equation (3.9) and then equation (B2) into the resulting expression yields 


x . (x) = - / h(x,y ) / k~ 1 (v,v' ;x)u! (v' ;x)dv'dv 

J V V v' 3 3 

= - / [ / h(x,y )k —1 (v ,v ' ;x)dv]u! (v ' ;x)dv ' 

v 1 v ^ ^ 

y B 

= - / sec\|) (y)[ f h(x,y )k— 1 (v ,v ' ;x)dv]u (x,y)dy , (B19) 

v b ' v j y nj 
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where, in the first line the original integral over y in equation (3.9) has been 
transformed into an integral over v, and in the last line the integral over v' in the 
second line has been transformed into an integral over y and equations (B5) and (B8) 
have been used. Substituting equation (3.15) into equation (B19) yields 

y B 

t (x) = -r / seo|i b (y) { / cos[^ b (y )-X (x,y )]kr 1 (v,v';x)dv} 
y A 

* u n j(x,y)dy 


y B 


+ / 


sea|i b (y)[ / (r Q -y v kos^y^kr 1 (v ,v^;x)dv]u n j (x,y)dy 


y R 

" ( t ) 

= - r Q / sea|» b (y)K^j (x,y)cos[i|i b (y)-A ^ (x.y^u^ (x,y)dy 
y A 


+ 


yB 

/ 

y A 


seci|/ b (y)[ / (r^y^cosi^y^kT 1 (v,Vy,*x)dv]u n j (x,y)dy , 


(B20) 


where the definition (B9) has been used in going to the second expression. The first 
lines of equations (3.16) and (B20) are identical. Thus, comparing the second lines 
of these two equations, we must have 


(r Q -y)K^ (x,y) = sect^Cy) / (r^y^cosi^y^kT 1 (v,v^;x)dv 


(B21 ) 


or 



(x,y) 


= { [ 1— (y/r )]cosi|> (y)} _ 1 

O D 


/ [l-(y /r )]cost(i (y )k— 1 (v ,v ' ;x)dv , 
vo b v j y 

v 


(B22) 
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in order that equations (3.16) and (B20) represent identical expressions for 
arbitrary deformations u n j(x,y). The integral with respect to v in equation (B22) is 
taken over the full length of line of contact of tooth pair j at zone of contact 
position x. 


To derive the integral equation whose solution yields K^j^(x,y), we apply the 
symmetry of equation (BIO) to equation (B22) which gives 


[ l-(y/r o )]costJ; b (y)K^' ) (x,y) = / kT 1 (v’.v ;x)[ l-^/r^] cosily )dv , (B23) 

J v J v 

which is of the same general form as equation (B2)* It follows from equations (Bl), 
(B2)» and (B23) that the local tooth-pair stiffness ^ (x ,y) can be obtained from 

the solution Pj^ V y> x ) ” — t * xe integral equation (Bl)*using for its left- 

hand side u f f (v;x) = u f ,(v»x) where 

j yj 


u 1 (v;x) ^ [ l-(y /r ) 1 cosij), (y ) 
yj L 'v o J r b 


(B24) 


and 


P uj (V y ;X) ~ [l-(y/r Q )]cosiJ) b (y)K^ , j ) (x,y) 


nj 


(B25) 


from which one obtains K^j^(x,y) _b^ 

K^(x,y) = [ l-(y/r o )]- 1 seo|) b (y)p^j(v y , ;x) , 


(B26) 


where the function * v ’ (y) is determined by equation (B6), 
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APPENDIX C 


DERIVATION OF EXPRESSION FOR h(x,y) GIVEN BY EQUATION (3.15) 

From Fig. 19 and the fact that the sum of the interior angles of a triangle 
is it radians , we have 


* = * ~ (f ' + b ) " X j = J + *b " X j * (C1) 

Applying the law of sines to the same triangle in Fig. 19 and using equation (Cl), we 
find 


r 




However , 


r o _ y 


sin(J + * b - 


(C2) 


sin (f " + b ) = cos \ 


(C3) 


and 


sin(y + i|» b - Xj) = cos(iJ> b -Xj) . 


(C4) 


Combining equations (C3) and (C4) with equation (C2) gives 

COS^k 

V 7 ’ ■ ( V y) 


(C5) 


or 


costly 

y ' ' r o ' co.(» t -X 3 ) ( V y) • 

Furthermore, from Fig. 19 and equation (Cl), we have 
JT = sin(ir-y) ■ sin[j - (ip b -Xj)] 

= cos(i|; b -Xj) . 


(C6) 


(C7) 
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Combining equations (C6) and (C7) yields the desired relationship for h: 


APPENDIX D 


INVERSION OF THE MESH STIFFNESS MATRIX K(x) 


From equations (2.10) and (3.37) through (3.42), it follows that the coefficient 
matrix K(x) in equation (3.40) can be written as 

(Dl) 

when the terms in the matrix (Dl) are appropriately identifed. Expressing the 
inverse C(x) of K(x) in terms of the cofactors (ref. 24, p. 17) of the matrix (Dl), 
it can be shown in a straightforward manner that the inverse of the matrix (Dl) can 
be expressed as 


K(x) = 


a l 

b. 


0 

0 


K -1 (x) = C(x) 




b 2 c 3 
” b l C 3 
_b 2 C l 


-a 2 e 3 


a l C 3 


a 2 C l 


0 

0 


a l b 2 -a 2 b l 


where the determinant K of the matrix (Dl) is 


K c^ (ajb 2 a 2 bj ) . 


(D2) 


(D3) 
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APPENDIX E 


DERIVATION OF EXPRESSION FOR MESH FORCE COMPONENTS IN GEAR COORDINATE SYSTEM 


We require the transformation 


{*<•>} = 


(El) 


that gives the three components x^*\ 1 = 1,2,3 of a vector originally expressed in 
terms of the components z^*\ ^ = 1*2,3 where the x ^ and z^ 1 ^ coordinate systems 
are illustrated in Fig. 12 or its mirror image, Fig. 13. The inverse of this 
transformation normally would be expressed as 


(•) 


{*'"} - [* (,) ] V°) , 

where [l' ^ ] 1 is the inverse of the matrix [i^* ^1 . However, since both the 
X^* ^ and z coordinate systems are orthogonal, we have (ref. 18, p. 103) 

[« <o r - [i <o ] . 


(E2) 


(E3) 


where is the transpose of the matrix . From equations (El) through (E3) 

and equation (A5), it follows that 


r (• ), 
[* 1 


cv 


(•) 


(•) (•) 
sy sv 

D 


(•) (•) 
qr b sv 




(•) 


"SY, 


(•) 


-sv 


(•) 


(•) (•) 


sy b cv 


(•) (•) 
CYb cv 


(E4) 


where c and s are abbreviations for cosine and sine operations as in equation (4.5). 
The angles v^*^ and Y^*^ are Illustrated in Figs. 12 and 13. Using equation (E4), 
one can readily verify that [2^*^] ■ [ i] , where [ i] is the identity matrix. 


Let W 


(•) 


z 


denote the six-component generalized force column matrix or vector 


described in the z. , i ■ 1,2,3 coordinate system analogous to the generalized force 
vector Q v * defined by equation (4.3) that describes the three torque components and 
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(•) 


three force components applied to gear (•) by its mating gear in the x ' ,i = 1,2,3 
coordinate system. Then can be obtained from W ^ by the transformation. 




(E5) 


where the matrix Gj^*^ is obtained from equation (E4): 


cv 


(*) . 


cy 


(•) 


) 

0 

- S v° > 

0 

0 


0 

»<*> 

■*b° 

sr<‘> cv<’> 

0 

0 


0 

sv ( -> 


cJ-J «<*> 

0 

0 


0 


0 

0 

CV (,) 

0 

-sv 

(•) 


0 

0 

sy<>> sv ( ° 

*b -) 


cv ( *> 


0 

0 

cy<*> sv ( *> 



CV (,) 


(E6) 


The six-component generalized force vector W ^ can be obtained from the three- 
component generalized mesh force vector W defined by equation (4.1) by a further 
matrix transformation. 


w<->=g 2 <->w. 


(E7) 


which, 

Q 


when combined with equation (E5), yields 


where 


«<•> c, ( - V-> 


(E8) 


(E9) 


We now evaluate the matrices \ G 2 ^*\ and G^ ^ for the various gearing con- 

figurations listed in Table 1 . 
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Driver clockwise with left-hand spiral or follower counterclockwise with left- 
hand spiral. Gear (1) shown in Fig. 3 is the clockwise rotating driver of the 
meshing gear pair shown there. According to the upper portion of the figure, gear 
(1) has a left-hand spiral (ref. 19, p. 241). The components of the total resultant 
force exerted on gear (1) by its mate, gear (2), are shown in Fig. 4. These 
components W t and W r and the resultant torque t taken about the zone of contact 
center jO^ are shown in Fig. 12. The force components W t and W r shown in Fig. 12 lie in 

the plane of contact, which is the plane occupied by the z^ , z 3 and y 1 , y 3 axes shown 

there. Torque x about the zone of contact center _0_ is represented as a vector normal 
to the plane of contact passing through 0_ with its direction governed by the right- 
hand rule, which is consistent with the sign convention described in the main text 
since the relationship between the positive direction of each p.^*^ and the 
corresponding x^ * , 1 * 1*2,3 in Fig. 12 is governed by the right-hand rule. 

To transform this mesh force system to the gear coordinate system x^*^ and 

i = 1,2,3 illustrated in Fig. 12, we must first express about the origin of 

the gear coordinate system located at the apex of the base cone shown in Fig. 12 the 

torque resulting from the total resultant mesh force. This translated torque 

vector lies along the z 2 axis, and has a value of x + r c W t * positive in the negative 

z 2 direction as shown in Fig. 12. Since W t and W f lie along the z ^ and z 3 axes, 

respectively, positive in the positive z directions, it follows for this case that 

the vector W \ defined above, is 
z 


-<T+r 0 W t )' 


w 

z ' 


y 


w 


w 


(E10) 


From the definition (1) of W, and equations (E7) and (E10), it follows directly 
that for this case the matrix G 2 ^*^ is 
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0 


0 


0 




(•) 


-1 -r 0 

o 

0 0 0 

0 1 0 

0 0 0 


0 


0 


1 


(Ell) 


The validity of equation (Ell) can be established by multiplying equation 

(Ell), by W, equation (4.1), which yields W ^ as prescribed by equation (E7). 

z 

Finally, the matrix is obtained by multiplying G^** and G 2 ^ as indicated 

by equation (E9): 


0 

0 

0 

-c Yb (,) 

"vV' ) 

0 


V T b ( ’ > 

0 

0 

«<•> 


0 


ST b ( ' , cv <-) 

0 


c ¥b <-V ) 


(E12) 


By comparing equations (4.5) and (E12), one can see that 


(E13) 


when the upper signs in the third column in the right-hand side of equation (4.5) are 
used, which is consistent with the sign convention described in Table 1. According 
to equation (E8), equation (E13) verifies equation (4.4) for the case of a driving 
gear rotating clockwise with left-hand spiral teeth as illustrated in Fig. 3. 
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Let us turn now to the case of a follower gear rotating counterclockwise with 
left-hand spiral teeth. This situation is depicted by the follower gear (2) shown in 
Fig. 3 except that the spiral angle is opposite to that shown in the upper portion of 
the figure. For this configuration, the total resultant force exerted on the 
follower gear (2) by the driver gear (1) also is correctly depicted by Fig. 4, and 
the force system illustrated in Fig. 12 also correctly depicts the forces exerted on 
the follower gear (2) by the driver gear (1). From this fact, it follows that the 
results provided by equations (El 2) and (El 3) also apply to the case of a follower 
gear rotating counterclockwise as depicted in Fig. 3, but with a left-hand spiral 
which is the opposite to that illustrated in the figure. 

Follower counterclockwise with right-hand spiral or driver clockwise with right- 
hand spiral. Gear (2) shown in Fig. 3 is the counterclockwise rotating follower of 
the meshing gear pair shown there. According to the upper portion of the figure, 
gear (2) has a right-hand spiral. Each force component W t , W r , and x exerted on gear 
(2) by its mate, gear (1), in this case is directed opposite to that shown in Fig. 4, 
since that figure illustrates the forces exerted on gear (1) by gear (2). 

The configuration shown in Fig. 12 is applicable to the follower gear (2), as 
well as the driver gear (1). When the above-described mesh force components w r , 
and x acting on gear (2) are sketched on Fig. 6, one obtains the force system 
illustrated in Fig. 12 with the exception that the component W r is directed opposite 
to that shown in Fig. 12. In this case, consequently, the six-component generalized 
force vector W which expresses the three torque components and three force 

components acting on the gear under consideration in the z^, i = 1,2,3 coordinate 
system, is equal to that given by equation (E10) with the exception that the sign of 
W r is reversed: 

r 

0 

-(x+ rW„) 

► . (E14) 


o t 


W<*> 
z 


0 

W t 

0 

-w 
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From the definition (4.1) and equations (E7) and (E14), it then follows that for the 
present case the matrix G 2 ^*^ is 


(•) 


0 

-1 

0 

0 

0 

0 


-r 

c 

0 

1 

0 

0 


0 

0 

0 

0 

0 

-1 


(E15) 


which is the same as equation (Ell) except for the sign change in the lowermost 
element in the third column. 

The matrix is obtained by multiplying the matrices described by equations 

(E6) and (E15) as indicated by equation (E9). This multiplication yields 


,(•) . 


-cy 


sy 


(•) 

b 

(•) 


" r o CY b 


r s Yu 
o 'b 


(•) 

(•) 


(•) 

(• ) 

cv 

sv 

(•) (•) 

(•) (•) 

sy b sv 

-sy k cv 

(•) (•) 

(•) (•) 

cy sv 

D 

-cy ^ cv 


(E16) 


0 


By comparing equations (4.5) and (E16), one can see that equation (E13) is valid in 
this case also when the lower signs in the third column in the right-hand side of 
equation (4.5) are used, which is consistent with the sign convention described in 
Table 1. 


We turn now to the case of a driver gear rotating clockwise but with right-hand 
spiral teeth. This case is depicted by the driver, gear (1), shown in Fig. 3 except 
that the spiral angle is opposite to that shown in the upper portion of the figure. 
For this case also, each force component W t , W r » and x exerted on the driver, gear 
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[ (1), by the follower, gear (2), is directed opposite to that shown in Fig. 4. Hence, 

■ for this case also, the system of forces exerted on the driver, gear (1), by the 
follower, gear (2), is represented by situation depicted in Fig. 12, with the 
exception that the component W r is directed opposite to that shown in Fig. 12. It 
follows that equations (E14), (E15), and (E16) are valid for this case also, which 
completes our proof of the validity of equations (4.4) and (4.5) for the last of the 
four gear element configurations enumerated in the first two lines of Table 1. 

Driver counterclockwise with right-hand spiral or follower clockwise with right- 
hand spiral. Each of the remaining four gear element configuration possibilities 
delineated in the last two lines in Table 1 represents an exact mirror image of one 
of the four configurations in the first two lines of the table. Specifically, the 
third and fourth gear element configurations delineated in each of the two columns in 
Table 1 represent, respectively, the exact mirror images of the first and second gear 
element configurations in the same column. 

The configuration shown in Fig. 7 is the mirror image of that shown in Fig. 3. 
The case of a counterclockwise rotating driver with right-hand spiral teeth is 
illustrated by Gear (1) of Fig. 7. The case of a clockwise rotating follower with 
right-hand spiral teeth is represented by gear (2) of Fig. 7, except that in this 
latter case the spiral angle is opposite to that shown in the upper portion of the 
figure. The forces exerted on the gear under consideration by its mating gear in 
each of these two cases are represented by the force system shown in Fig. 4. For 
both cases, this force system also is shown in Fig. 13. 

Figure 13 is an exact mirror image of Fig. 12. From Fig. 13, it follows that 

the expression for the generalized force vector W ^ also is given by equation (E10) 

z 

for the two gear element configurations under present consideration. Hence, the 

(•) (•) 

expressions given by equations (Ell) and (E12) for the matrices G^ and G , 
respectively, also are applicable to the present two cases, from which it follows 
that equations (E13), (4.4), and (4.5) also are verified when the upper signs in the 
third column in the right-hand side of equation (4.5) are used, which is consistent 
with the sign convention described in Table 1. However, it is important to recognize 
that for the present two cases, and the two to be described below, the sign 
convention for each i = 1,2,3 describing the relationship between each rotation 
component p^* \ and hence, and the displacement component x^*^ is governed 

by the left-hand rule, as may be seen in Fig. 8 or Fig. 3. 
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Follower clockwise with left-hand spiral or driver counterclockwise with left- 
hand spiral. The first of these two cases is represented by gear (2) of Fig. 7, and 
the second is represented by gear (1) of the same figure except that in this latter 
case the spiral angle is opposite to that shown in the upper portion of the figure. 

In each of these two cases , each of the force components exerted on the gear under 
consideration by its mating gear is directed opposite to that illustrated in Fig. 4. 
These force components exerted on the gear under consideration are shown in Fig. 13 
except that, for the present two cases, the direction of W r is opposite to that shown 
in Fig. 13. It follows that, for the present two cases, the generalized force 
vector W^ # ^ is given by equation (E14); hence, the matrices G^*^ and G^ # ^ are given 
by equations (E15) and (E16), respectively. Therefore, equation (E13) is valid for 
these two cases when the lower signs in the third column in the right-hand side of 
equation (4,5) are used, which is consistent with the sign convention described in 
Table 1. 
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APPENDIX F 


SIGN OF AjU.y^) 

The definition (5.54) expresses Aj(x,y e j) as 


. 4 Vij ' r o b i < *y-»u ) 

j a b - b a 

12 12 


ab - ba = a b fl •) 

12 12 1 2 a b_ J 


Combining equations (Fl) and (F2) gives 


. L ViVUVij -'.YVV' 

Vlj 


b a \-l 
1 2 \ 1 


- 1 - 


, S. . ^ i • 

b 2 lj o a j 2j lj 


i From equation (F3) it is clear that if bj = 0 and (gn/b 2 ) > 0, then Aj > 0. By 

<v ( t ) 

| comparing equation (5.3) with equation (3.38), we see that b = K (x;* ) , which 

2 c— 1 b 

from equations (3.36a) and (3.32a) may be seen always to be positive. Furthermore, 
r it follows from equation (5.22) that g^j always is positive; hence, (gjj/b 2 ) > 0. 
Therefore, it is necessary to show only that we can set b^ = 0. 

The coordinate system we have used was devised to keep equations (3.37) and 
(3.38) relatively uncoupled. In fact, if bj ■ 0, equation (3.38), can be solved for 
independently of the value of y. From equations (5.3) and (3.38), it follows that 

1 b l = (x;»|) ,X ) + r K^? (x;i|> ,X )] , (F4) 

yc x c b j o c 1(1 - c ) b j 

I 

where the notation is described by equations (3.32g), (3.32d), (3.36g), and (3.36d). 

: ~(t) 

ml i rr v / / . . \ . . . /«,« \ * _ j ^ j j _ J 1 a. _ n 


The term K 


c-l(l-c)'” ,r b’ j 


(x;i|> ,X ) in equation (F4) is positive and independent of r Q . 


Therefore, the magnitude of the positive term r K . (x;t|» ,X ) can be adjusted by 

o c-i(l-c) b j 
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~(t ) 

changing the value of r . The term K (x;i|> ) can be either positive or 

u yc i c b j 

negative depending on the value of r Q . However, careful examination of equations 

(3.32g), (3.32d), and Fig. 11 shows that r Q always can be chosen so that 

~(t) yc 1 c 

(x;ij> ,A ) is made negative with magnitude equal to r K (x;ij> ,X ), which will 

D J o c 1 (l-c) b j 

set b| = 0. Increasing r Q will decrease the value of bj, and decreasing r Q will 

increase the value of bj. 

The value of r Q controls the position of the origin of the coordinate y, as one 
can see from Fig. 11. We have chosen the value of r Q to position the origin of the 
coordinate y in the middle of the nominal contact range F, as shown in Fig. 11. 

This choice tends to keep the magnitude of bj very small, as careful examination of 
equations (3.32g), (3.32d), (3.36g), (3.36d), (F4), and Fig. 11 will show. Thus, 
with the choice of r Q shown in Fig. 11, we can expect Aj(x,y e j) usually to be 
positive, requiring in these cases no readjustment in the value of r Q . 
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